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1 Introduction 

Understanding the fundamental structure of matter is one of the central goals 
of scientific research. In the closing decades of the twentieth century, physicists 
developed a beautiful theory, Quantum Chromodynamics (QCD), which explains 
all of strongly interacting matter in terms of point-like quarks interacting by the 
exchange of gauge bosons, known as gluons. The gluons of QCD, unlike the 
photons of QED, can interact with each other. The color force which governs 
the interaction of quarks and gluons is responsible for more than 99% of the 
observable mass in the physical universe and explains the structure of nucleons 
and their composite structures, atomic nuclei, as well as astrophysical objects 
such as neutron stars. 

During the last 30 years, experiments have verified QCD quantitatively in 
collisions involving a very large momentum exchange between the participants. 
These collisions occur over very short distances much smaller than the size of the 
proton. In these experiments, the confined quarks and gluons act as if they are 
nearly free pointlike particles and exhibit many properties that arc predicted by 
perturbativc QCD (pQCD). This experimental phenomenon was first discovered 
in deeply inelastic scattering (DIS) experiments of electrons off nucleons. The 
discovery resulted in the 1990 Nobel Prize in Physics being awarded to Friedman, 
Kendall and Taylor. The phenomenon, that quarks and gluons are quasi-free at 
short distances, follows from a fundamental property of QCD known as asymp- 
totic freedom. Gross, Politzer and Wilczek, who first identified and understood 
this unique characteristic of QCD were awarded the 2004 Nobel Prize in Physics. 

When the interaction distance between the quarks and gluons becomes compa- 
rable to or larger than the typical size of hadrons, the fundamental constituents 
of the nucleon are no longer free. They are confined by the strong force that 
does not allow for the observation of any "colored" object. In this strong cou- 
pling QCD regime, where most hadronic matter exists, the symmetries of the 
underlying quark-gluon theory are hidden, and QCD computations in terms of 
the dynamical properties of quarks and gluons are difficult. A major effort is 
underway worldwide to carry out ab initio QCD calculations in the strong QCD 
regime using Monte-Carlo simulations on large scale computers. 

The experimental underpinnings for QCD are derived from decades of work 
at the CERN, DESY, Fcrmilab and SLAC accelerator facilities. Some highlights 
include the determination of the nucleon quark momentum and spin distributions 
and the nucleon gluon momentum distribution, the verification of the QCD pre- 
diction for the running of the strong coupling constant a^, the discovery of jets, 
and the discovery that quark and gluon momentum distributions in a nucleus 
differ from those in a free nucleon. 

However, thirty years after QCD has been established as the Standard Model of 
the strong force, and despite impressive progress made in the intervening decades, 
understanding how QCD works in detail remains one of the outstanding issues in 
physics. Some crucial open questions that need to be addressed are listed below. 
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- What is the gluon momentum distribution in the atomic nucleus? 

QCD tells us that the nucleon is primarily made up of specks of matter 
(quarks) bound by tremendously powerful gluon fields. Thus atomic nuclei 
are primarily composed of glue. Very little is known about the gluon mo- 
mentum distribution in a nucleus. Determining these gluon distributions 
is therefore a fundamental measurement of high priority. This quantity is 
also essential for an understanding of other important questions in hadronic 
physics. For example, the interpretation of experiments searching for a de- 
confined quark-gluon state in relativistic heavy ion collisions is dependent 
on the knowledge of the initial quark and gluon configuration in a heavy 
nucleus. This will be especially true for heavy ion experiments at the Large 
Hadron Collider (LHC) at CERN. Further, there are predictions that glu- 
onic matter at high parton densities has novel properties that can be probed 
in hard scattering experiments on nuclear targets. Hints of the existence of 
this state may have been seen in Deuteron-Gold experiments at the Rela- 
tivistic Heavy Ion Collider (RHIC) at Brookhaven. 

- How is the spin structure of the nucleon understood to arise from the quark 

and gluon constituents ? 

High energy spin-dependent lepton scattering experiments from polarized 
nucleon targets have produced surprising results. The spins of the quarks 
account for only about 20% of the spin of the proton. The contribution of 
the gluons may be large. Dramatic effects arc predicted for measurements 
beyond the capability of any existing accelerator. There are hints from 
other experiments that the contribution of orbital angular momentum may 
be large. 

- Testing QCD 

It is imperative to continue to subject QCD to stringent tests because there 
is so much about the theory that remains a mystery. QCD can be tested 
in two ways: one is by precision measurements, and the other is by looking 
for novel physics which is sensitive to the confining properties of the theory. 
Both of these can be achieved at a high luminosity lepton-ion collider with 
a detector that has a wide rapidity and angular coverage. An example of 
precision physics is the Bjorken Sum Rule in spin-dependent Icpton scatter- 
ing from a polarized nucleon. This fundamental sum rule relates inclusive 
spin-dependent lepton scattering to the ratio of axial to vector coupling con- 
stants in neutron /?-decay. Present experiments test it to about ±10%: it 
would be highly desirable to push these tests to about ±1%. Further, with 
lattice QCD expected to make substantial progress in the ability to make 
ah initio QCD calculations during the next decade, precise measurements 
of the calculable observables will be required. An example of a physics 
measurement sensitive to confinement is hard diffraction, where large mass 
final states are formed with large "color-less" gaps in rapidity separating 
them from the hadron or nucleus. At the Hadron Electron Ring Accelerator 
(HERA) at DESY, roughly 10% of events are of this nature. The origins 
of these rapidity gaps, which must be intimately related to the confining 
properties of the theory, can be better understood with detectors that are 
able to provide detailed maps of the structure of events in DIS. 

This article motivates and describes the next generation accelerator required 
by nuclear and particle physicists to study and test QCD, namely a polarized 
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lepton-ion collider. The basic characteristics of the collider are motivated as 
follows: 

- lepton beam 

The lepton probe employs the best understood interaction in nature (QED) 
to study hadron structure. Electrons and positrons couple directly to the 
quarks. The experimental conditions which maximize sensitivity to valence 
and sea quarks as well as probe gluons are well understood. Further, the 

availability of both positron and electron beams will enable experiments 
that arc sensitive to the exchange of the parity violating Z and W-bosons. 

- range of center- of-mass (CM) energies 

To cleanly interact with quarks, a minimum center-of-mass (CM) energy of 
about 10 GeV is required. To explore and utilize the powerful evolu- 
tion equations of QCD, CM energies of order 100 GeV are desirable. This 
consideration strongly motivates the collider geometry. 

- high luminosity 

The QED interaction between the lepton probe and the hadron target is 
relatively weak. Thus precise and definitive measurements demand a high 
collision luminosity of order 10^^ nucleons cm~^ s^^. 

- polarized beams 

Polarized lepton and miclcon beams are essential to address the central 
question of the spin structure of the nucleon. Both polarized proton and 
neutron (effectively polarized or ^He) are required for tests of the fun- 
damental Bjorken Sum Rule. The polarization direction of at least one of 
the beams must be reversible on a rapid timescale to minimize systematic 
uncertainties. 

- nuclear beams 

Light nuclear targets are useful for probing the spin and flavor content of 
parton distributions. Heavy nuclei are essential for experiments probing the 
behavior of quarks and gluons in the nuclear medium. 

- detector considerations 

The collider geometry has a significant advantage over fixed-target exper- 
iments at high energy because it makes feasible the detection of complete 
final-states. A central collider detector with momentum and energy mea- 
surements and particle identification for both leptons and hadrons will be 
essential for many experiments. Special purpose detectors that provide wide 
angular and rapidity coverage will be essential for several specific measure- 
ments. 

These considerations constrain the design parameters of the collider to be a 
5 to 10 GeV energy electron (or positron) beam colliding with a nucleon beam 
of energy 25 GeV to 250 GeV. The collider is anticipated to deliver nuclear 
beams of energies ranging from 20 — 100 GeV/nucleon. The lepton and nucleon 
beams must be highly polarized and the collision luminosity must be of order 
10^^ nucleons cm~^ s~^. The proposed eRHIC design (described in section 4) 
realizes the required specifications in a cost effective and timely way by using the 
existing RHIC facility at BNL. The characteristics of eRHIC are well beyond the 
capability of any existing accelerator, as is clear from Fig. 1. 

By delivering high energies to the collision, the collider provides an increased 
range for investigating quarks and gluons with small momentum fraction (x) and 
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Figure 1: The center-of-mass energy vs. luminosity of the proposed Electron-Ion 
Collider eRHIC compared to other lepton scattering facilities. 



for studying their behavior over a wide range of momentum transfers (Q^). In 
deeply inelastic scattering, the accessible values of the Bjorken variable x (defined 
in section 2) are limited by the available CM energy. For example, collisions 
between a 10 GeV lepton beam and nuclear beams of 100 GeV/nucleon provide 
access to values of x as small as 3 x 10"^ for Q2 ^ 1 GeV^. In a fixed-target 
configuration, a 2.1 TeV lepton beam would be required to produce the same 
CM energy. Figure 2 shows the x-Q^ range possible with the proposed eRHIC 
machine and compares that range to the currently explored kinematic region. 

In this article, the scientific case and accelerator design for a new facility to 
study the fundamental quark and gluon structure of strongly interacting matter 
is presented. Section 2 describes the current understanding of the quark and 
gluon structure of hadrons and nuclei. Section 3 presents highlights of the sci- 
entific opportunities available with a lepton-ion collider. Section 4 describes the 
accelerator design effort and section 5 describes the interaction region and eRHIC 
detector design. 

2 Status of the Exploration of the Partonic Structure of Hadrons 
and Nuclei 

This section will summarize our current understanding of the partonic structure 
of hadrons and nuclei in QCD, accumulated during the past three decades from 
a variety of deeply inelastic and hadronic scattering experiments. We will also 
comment on what new information may become available from DIS as well as from 
RHIC and other experimental facilities around the world before a future electron- 
ion collider starts taking data. We will outline the status of our knowledge on 
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Figure 2: The x-Q^ range of the proposed lepton-ion cohider at Brookhaven 
National Laboratory (eRHIC) in comparison with the past & present experimen- 
tal DIS facilities. The left plot is for polarized DIS experiments, and the right 
corresponds to the unpolarized DIS experiments. 



i) the parton distributions in nucleons, ii) spin and flavor distributions in the 
nucleon, iii) nuclear modifications to the inclusive nucleon distributions such as 
the European Muon Collaboration (EMC) effect and quark and gluon shadowing, 
iv) color coherent phenomena in nuclei that probe the space-time structure of 
QCD such as color transparency and opacity, partonic energy loss and the px 
broadening of partons in media. In each case, we will outline the most important 
remaining questions and challenges. These will be addressed further in Section 3. 

2.1 Deeply-inelastic scattering 




Figure 3: Deeply- inelastic lepton- nucleon scattering mediated by virtual photon 
exchange. 



The cross-section for the inclusive deeply inelastic scattering (DIS) process 
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shown in Fig. 3 can be written as a product of the leptonic tensor C^i, and the 
hadronic tensor W'^'^ as 

^ cx j(:^,{k,q,s)W^'''{P,q,S), (1) 

where one defines the Lorentz invariant scalars, the famous Bjorken variable 
X = —q^/2Pq, and y = Pq/Pk. Note that as illustrated in Fig. 3, k {k') is the 
4-momentum of the incoming (outgoing) electron, P is the 4-momentum of the 
incoming hadron, and q = k — k' is the 4-momentum of the virtual photon. The 
center of mass energy squared is s = (P + A;)^. From these invariants, one can 
deduce simply that xy ^ Q'^/s, where = —q^ > 0. 
The hadronic tensor can be written in full generality as 

W''''{P,q,S) = ^J d^ze"i'{P,S\ [Mz),J.m \P.S) = -g''''F^{x,Q^) 



F2{x,Q')-iet'-P''^F^{x,Q') + ie^'^'^^qp 



Pq ' - ' 2Pq ' - ' ^'^ [Pq 



2\ 



S.{P q) - P.{S q) 2^ 

+ lp^2 92{X,Q ) 



+ 



p,S- + S^^P- Sq , 
2Pq {Pqf 



g3{x,Q ) 



+ (p^^'^^' Q') - -f-^g^- g^{x, Q') . (2) 

The Fi are referred to as the "unpolarized" structure functions, whereas the 
gi are the "spin-dependent" ones, because their associated tensors depend on 
the nuclcon spin vector S^. Note that parity- violating interactions mediated by 
electroweak boson exchange are required for F^, g^, g^, g^ to contribute. 

Inserting Eq. 2 and the straightforwardly calculated leptonic tensor into Eq. 1, 
one obtains the DIS cross section in terms of the structure functions. If one 
averages over the hadronic spins and restricts oneself to parity conserving (for 
^ M'^) electron-nucleon scattering alone, one finds the simple expression 



da 27ragjjj 



(l + (1 - yf) F2{x, Q2) - y-^FLix, Q2)] . (3) 



Here, Oem is the coupling constant of Quantum Electrodynamics and F^ is the 
"longitudinal" structure function, defined by the relation Fl = F2 — 2xFi. 

In the leading logarithmic approximation of QCD the measured structure func- 
tion F2{x,Q'^) can be written as 

F2{x, Q2) = (^«(^' + ^^(a^' Q')) ' (4) 

q=u,d,s,c,b,t 



where q{x,Q ) {q{x,Q )) is the probability density for finding a quark (anti- 

°q 



quark) with momentum fraction a; at a momentum resolution scale Q^; Cg is the 



quark charge. 

In the simple parton model one has Bjorken scaling, F2{x,Q^) — F2{x). The 
"scaling violations" seen in the Q^-dependence of F2{x,Q'^) arise from the fact 
that QCD is not a scale invariant theory and has an intrinsic scale Aqcd ~ 200 
MeV. They are only logarithmic in the Bjorken limit of — ^ 00 and s — ^ 00 
with X ~ Q^/s fixed. As one moves away from the asymptotic regime, the scal- 
ing violations become significant. They can be quantitatively computed in QCD 
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Figure 4: The plot on the left shows the world data on F2 as a function of for 
fixed values of x. On the right we show the converse: F2 as a function of x for 
fixed values of Q^. From (8). 



perturbation theory using for example the machinery of the operator product 
expansion and the renormalization group. The result is most conveniently sum- 
marized by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
equations for the parton densities (1,2): 

:7r^ (')ix, Q') = { p^^h' i p^^S"- "I ) ^ f ) (x, Q^) , (5) 

dlnQH^y \Pgq{as,x) Pgg{as,x) J \gJ\'^J' 

where (8) denotes a convolution, and the Pij are known as "splitting functions" (2) 
and are evaluated in QCD perturbation theory. They are now known to three-loop 
accuracy (3). The evolution of the quark densities q, q involves the gluon density 
g{x,Q'^). The physical picture behind evolution is the fact that the virtuality 
of the probe sets a resolution scale for the partons, so that a change in 
corresponds to a change in the parton state seen. The strategy is then to 
parameterize the parton distributions at some initial scale = Qq, and to 
determine the parameters by evolving the parton densities to (usually, larger) 
and by comparing to experimental data for F2{x, Q^). 

The pioneering DIS experiments, which first measured Bjorken scaling of F2, 
were performed at SLAC (4). However, because of the (relatively) small energies, 
these experiments were limited to the region of x > 0.1. With the intense muon 
beams of CERN and Fermilab, with energies in excess of 100 GeV, the DIS 
cross-section of the proton was measured down to and below x ~ 10^^ (5). In 
the 1990's, the HERA collider at DESY extended the DIS cross-section of the 
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proton to below x = 10"'* (6, 7). The current experimental determination of 
^proton Q^) extends over 4 orders of magnitude in x and Q^. This is shown in 
Fig. 4. The left panel in Fig. 4 shows next-to-leading order (NLO) QCD global 
fits by the ZEUS and HI detector collaborations at HERA to F2 as a function of 
for the world DIS data. The data and the QCD fit are in excellent agreement 
over a wide range in x and Q^. In the right panel of Fig. 4, the x dependence 
of F2 is shown for different bins in Q^. The rapid rise in F2 with decreasing x 
reflects the sizeable contribution from the sea quark distribution at small x. 




Figure 5: The valence (up and down) quark, sea quark and gluon distributions 
plotted as a function of x for fixed = 10 GeV^. Note that the sea and glue 
distributions are scaled down by a factor of 1/20. From (8). 



In Fig. 5 we show the valence up and down quark distributions as well as the 
gluon and sea quark distributions extracted by the HI and ZEUS collaborations 
as functions of x for fixed = 10 GeV^. The valence parton distributions 
are mainly distributed at large x whereas the glue and sea quark distributions 
dominate hugely at small x. Indeed, the gluon and sea quark distributions are 
divided by a factor of 20 to ensure they can be shown on the scale of the plot. 
Already at x ~ 0.1, the gluon distribution is nearly a factor of two greater than 
the sum of the up and down quark valence distributions. 

As follows from Eq. 5, the gluon distribution in DIS may be extracted from 
scaling violations of F2: xg{x, Q^) oc ^^qx^q-i ^ - As one goes to low Q^, xg{x, Q"^) 
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becomes small, and some analyses find a preference for a negative gluon distri- 
bution at low X, modulo statistical and systematic uncertainties (9,10). This is 
in principle not a problem in QCD beyond leading order. However, the resulting 
longitudinal structure function Fl also comes out close to zero or even negative 
for ~ 2 GeV^ ^, which is unphysical because Fl is a positive-definite quan- 
tity. A likely explanation for this finding is that contributions to Fl that are 
suppressed by inverse powers of are playing a significant role at these values 
of (11)- These contributions are commonly referred to as higher twist effects. 




Figure 6: The plot on the left shows the longitudinal structure function Fl as a 
function of x for different bins (9). On the right (from Ref. (12)) is a plot of 
the virtual photon-proton cross-section plotted as a function of r = Q'^ /Q^- See 
text for further explanation. 



It has been shown recently (12) that the HERA data on the virtual photon- 
proton cross-section {a"'*'P = Att'^ aemF2{x , Q'^) / Q'^) , for all x < 10^^ and 0.045 < 
< 450 GeV^, exhibit the phenomenon of "geometrical scaling" shown in 
Fig. 6. The data are shown to scale as a function of r = Q'^/QI, where Q^{x) = 
Qlixo/x)-^ with Ql = 1 GeV2, xq = 3-10-^ and A « 0.3. The scale is called 
the saturation scale. Geometrical scaling, although very general, is realized in a 
simple model, the Golec-Biernat-Wiisthoff model which includes all twist contri- 
butions (13). The model (and variants) provide a phenomenological description 
of the HERA data on diffractive cross-sections and inclusive vector meson pro- 
duction (14-18). The saturation scale and geometrical scaling will be discussed 
further in Section 3. 

^The leading twist expression for Fl is simply related to as xg{x,Q'^). 
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2.2 Spin structure of the nucleon 

2.2.1 What we have learned from polarized DIS 

Spin physics has played a prominent role in QCD for several decades. The field 
has been driven by the successful experimental program of polarized deeply- 
inelastic lepton- nucleon scattering at SLAC, CERN, DESY and the Jefferson 
Laboratory (19). A main focus has been on measurements with longitudinally 
polarized lepton beam and target. For leptons with helicity A scattering off 
nucleons polarized parallel or antiparallel to the lepton direction, one has (20) 

^fedQ2 d^dQ2~ ^ C{G^,Ga,X) [Xxy{2-y)gi + {l-y)gi + xy'^ g5 

(6) 

where C{Gv, Ga, A) are factors depending on the vector and axial couplings of the 
lepton to the exchanged gauge boson. The terms involving 34 and (75 in Eq. 6 are 
associated with Z and W exchange in the DIS process and violate parity. In the 
fixed-target regime, pure-photon exchange strongly dominates, and scattering off 
a longitudinally polarized target determines gi. Figure 7 (left) shows a recent 
compilation (21) of the world data on gi{x, Q"^), for proton, deuteron, and neutron 
targets. Roughly speaking, gi is known about as well now as the unpolarized F2 
was in the mid-eighties, prior to HERA. Figure 7 (right) shows the measured 
Q^-dependence of gw the predicted scaling violations are visible in the data. 
In leading order of QCD, gi can be written as 
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Figure 7: Left: world data on the spin structure function gi as compiled and 
shown in (21). Right: gi{x, Q^) as a function of for various x. The curves are 
from a phenomenological fit. Taken from (22). 
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where 

^q = qZ - {q = u,d,s,...) , (8) 

(q^) denoting the number density of quarks of same (opposite) hehcity as the 
nucleon. Clearly, the Ag(x, Q^), Ag(x, Q^) contain information on the nucleon 
spin structure. Also in the spin-dependent case, QCD predicts Q^-dependence of 
the densities. The associated evolution equations have the same form as Eq. 5, 
but with polarized splitting functions (2,23,24). Also, the spin-dependent gluon 
density Ag, defined in analogy with Equation 8, appears. 

The results of a recent QCD analysis (25) of the data for gi{x, Q^) in terms of 
the polarized parton densities are shown in Fig. 8. The shaded bands in the figure 
give estimates of how well we know the distributions so far. As can be seen, the 
valence densities are fairly well known and the sea quark densities to some lesser 
extent. This analysis (25) assumes flavor-SU(3) symmetry for the sea quarks; the 
actual uncertainties in the individual sea distributions are much larger. Finally, 
Fig. 8 also shows that we know very little about the polarized gluon density. A 
tendency toward a positive Ag is seen. It is not surprising that the uncertainty in 
Ag is still large: at LO, Ag enters only through the Q^-evolution of the structure 
function gi. Because all polarized DIS experiments thus far have been with fixed 
targets, the lever arm in has been limited. This is also seen in a comparison 
of Fig. 7 with Fig. 4. 




Figure 8: Recent analysis of polarized parton densities of the proton. Taken from 
(25) ("BB"). The additional curves represent the central fits from the analyses 
of (26) ("GRSV") and (27) ("AAC"). 



A particular focus in the analysis of gi has been on the integral Ti(Q^) = 
Jq gi{x,Q'^)dx. Ignoring QCD corrections, one has from Eq. 7: 

Fi = ^AAs + ^AAs + ^AS , (9) 
12 00 9 
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where 



AS 



AU + AU + AV + AV + AS + AS , 

AU+ AU - AV - AV , 

AU+ AU + AV + AV-2 (AS + AS) 



AAs 



(10) 



with AQ = Jq Aq{x, Q'^)dx, which does not evolve with at lowest order. The 
flavor non-singlet combinations A^, turn out to be proportional to the nucleon 
matrix elements of the quark non-singlet axial currents, (P, S \ q^^ 'j^ Xi q\ P, S) . 
Such currents typically occur in weak interactions, and by SU(3) rotations one 
may relate the matrix elements to the /?-decay parameters F, D of the baryon 
octet (28,29). One finds AA3 = F + D = gA = 1-267 and AAs = 3F - D 0.58. 
The first of these remarkable connections between hadronic and DIS physics 
corresponds to the famous Bjorken sum rule (28), 



where the superscripts p and n denote the proton and neutron respectively. The 
sum rule has been verified experimentally with about 10% accuracy (19). The 
QCD corrections indicated in Equation 11 are known (30) through ©(a^). Assum- 
ing the validity of the sum rule, it can be used for a rather precise determination 
of the strong coupling constant (31). 

Determining Fi from the polarized-DIS data, and using the information from 
/3-decays on A A3 and A^g as additional input, one may determine AE. This 
quantity is of particular importance because it measures twice the quark spin 
contribution to the proton spin. The analysis reveals a small value AS ~ 0.2. The 
experimental finding that the quarks carry only about 20% of the proton spin 
has been one of the most remarkable results in the exploration of the structure 
of the nucleon. Even though the identification of nucleon with parton helicity 
is not a prediction of QCD (perturbative or otherwise) the result came as a 
major surprise. It has sparked tremendous theoretical activity and has also been 
the motivation behind a number of dedicated experiments in QCD spin physics, 
aimed at further unraveling the nucleon spin. 

A small value for AS also implies a sizable negative strange quark polarization 
in the nucleon, AS + AS —0.12. It would be desirable to have independent ex- 
perimental information on this quantity, to eliminate the uncertainty in the value 
for AS due to SU(3) breaking effects in the determination of AAs from baryon 
j3 decays (32). More generally, considering Fig. 8, more information is needed 
on the polarized sea quark distribution functions and their flavor decomposition. 
Such knowledge is also very interesting for comparisons to model calculations of 
nucleon structure. For example, there have been a number of predictions (33) for 
the An — Ad. Progress toward achieving a full flavor separation of the nucleon 
sea has been made recently, through semi-inclusive measurements in DIS (SIDIS) 
(34,35). Inclusive DIS via photon exchange only gives access to the combinations 
Aq + Aq, as is evident from Equation 7. If one detects, however, a hadron in the 
flnal state, the spin-dependent structure function becomes 



Here, the D^{z) are fragmentation functions, with z = E^/u, where is the 
energy of the produced hadron and v the energy of the virtual photon in the Lab 



F? - F^ = i A^3 [1 + 0{as)\ = \gA [l + C(«,) 



(11) 



(12) 
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frame. Fig. 9 shows the latest results on the flavor separation by the HERMES 
collaboration at HERA (35). Uncertainties are still fairly large; unfortunately, 
no further improvements in statistics are expected from HERMES. The results 
are not inconsistent with the large negative polarization of Au = Ad = As in 
the sea that has been implemented in many determinations of polarized parton 
distributions from inclusive DIS data (see, e.g., the curves in Fig. 8). On the 
other hand, there is no evidence for a large negative strange quark polarization. 
The results have sparked much renewed theory activity on SIDIS (36). We note 
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Figure 9: Recent HERMES results (35) for the quark and antiquark polarizations 
extracted from semi-inclusive DIS. 



that at RHIC production will be used to determine Au, Au, Ad, Ad with 
good precision, exploiting the parity-violating couplings of the W to left-handed 
quarks and right-handed antiquarks (37). Comparisons of such data taken at 
much higher scales with those from SIDIS will be extremely interesting. 

A measurement of Fi obviously relies on an estimate of the contribution to 
the integral from x outside the measured region. The extrapolation to small 
X constitutes one main uncertainty in the value of AS. As can be seen from 
Fig. 7, there is not much information on gi{x,Q'^) sX x < 0.003. In addition, 
the data points at the smaller x also have values that are below the DIS 
regime, making it conceivable that the "higher-twist" contributions to gi{x,Q'^) 
are important and contaminate the extraction of AS. About half of the data 
points shown in Fig. 7 are from the region < 4 GeV^, W"^ = Q^(l — x)/x < 10 
GeV^, which in the unpolarized case is usually excluded in analyses of parton 
distribution functions. Clearly, measurements of polarized DIS and SIDIS at 
smaller x, as well as at presently available x, but higher Q^, will be vital for 
arriving at a definitive understanding of the polarized quark distributions, and 
of AS in particular. 
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2.2.2 Contributors to the nucleon spin 

The partons in the nucleon have to provide the nucleon spin. When formulating 
a "proton spin sum rule" one has in mind the expectation value of the angular 
momentum operator (38,39), 

i = (P, 1/2 I J3 I P, 1/2) = (P, 1/2 I y d3 [f X f ] J P, 1/2) , (13) 

where = with T the QCD energy-momentum tensor. Expressing the 
operator in terms of quark and gluon operators, one may write: 



^ = ^AS + AG(Q2) + L,(Q2) + LgiQ^) , (14) 



where AG(Q^) = Ag{x,Q^) is the gluon spin contribution and the Lg^g cor- 
respond to orbital angular momenta of quarks and gluons. Unlike AS, AG and 
Lq^g depend on the resolution scale already at lowest order in evolution. The 
small size of the quark spin contribution implies that wc must look elsewhere for 
the proton's spin: sizable contributions to the nucleon spin should come from 
AG and/or Lg^g. 

Several current experiments are dedicated to a direct determination of Ag{x, Q^). 
High-transverse momentum jet, hadron, and photon final states in polarized pp 
scattering at RHIC offer the best possibilities (37). For example, direct access to 
Ag is provided by the spin asymmetry for the reaction pp — jX, owing to the 
presence of the QCD Compton process qg — > -jq. The Spin Muon Collaboration 
(SMC) and COMPASS fixed-target experiments at CERN, and the HERMES 
experiment at DESY, access Ag{x,Q'^) in charm or high-pr hadron pair final 
states in photon-gluon fusion 7*3 — > qq (40) . Additional precision measurements 
with well established techniques will be needed to determine the integral of the 
polarized gluon distribution, particularly at lower x. 

Orbital effects are the other candidate for contributions to the proton spin. 
Close analysis of the x x T matrix elements in Eq. 13 revealed (38) that they 
can be measured from a wider class of parton distribution functions, the so- 
called generalized parton distributions (GPD) (41). These take the general form 
{p + A\Oq^g\p) , where Og^g are suitable quark and gluon operators and A is some 
momentum transfer. The latter is the reason that the GPDs are also referred to 
as "off-forward" distributions. The explicit factor x in Equation 13 forces one 
off the forward direction, simply because it requires a derivative with respect to 
momentum transfer. This is in analogy with the nuclcon's Pauli form factor. 
In fact, matrix elements of the above form interpolate between DIS structure 
functions and elastic form factors. 

To be more specific (42), the total (spin plus orbital) angular momentum con- 
tribution of a quark to the nucleon spin is given as (38) 

Jg = I lim [dxx \Hg{x,^,A^) + Eg{x,C,A^)\ . (15) 

Here, ^ = A"'"/P^, where the light-cone momentum A"*" = A'^ + A^, and likewise 
for P"*" . Hg , Eg are defined as form factors of the matrix element 
/ dye^y^ (P' I i^^{y) ip+{0) |P). Hg reduces to the ordinary (forward) quark 
distribution in the limit A — 0, Hq{x, 0, 0) = q{x), whereas the first moments (in 
x) of Hg and Eg give the quark's contributions to the nucleon Dirac and Pauli 
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form factors, respectively. In addition, Fourier transforms of Hq,Eq with respect 
to the transverse components of the momentum transfer A give information on 
the position space distributions of partons in the nucleon (43), for example: 

Hg{x, ^ = 0, -Ai) = J d%e-'^^ 'qix, b) . (16) 

q{x, b) is the probability density for finding a quark with momentum fraction x 
at transverse distance b from the center. It thus gives a transverse profile of the 
nucleon. GPDs, therefore, may give us remarkably deep new insight into the 
nucleon. 

The classic reaction for a measurement of the Hg, Eg is "deeply virtual Comp- 
ton scattering (DVCS)", 7*p (38). It is the theoretically best explored and 

understood reaction (44). Next-to- leading order calculations are available (45). 
The GPDs contribute to the reaction at amplitude level. The amplitude for 
DVCS interferes with that for the Bethe-Heitler process. The pure Bcthc-Heitler 
part of the differential ep cross section is calculable and can in principle be sub- 
tracted, provided it does not dominate too strongly. Such a subtraction has been 
performed in DVCS measurements at small x by HI and ZEUS (46). A different 
possibility to eliminate the Bethe-Heitler contribution is to take the difference 
of cross sections for opposite beam or target polarization. In both cases, con- 
tributions from Compton scattering and the Compton-Bethe-Heitler interference 
survive. The cleanest separation of these pieces can be achieved in experiments 
with lepton beams of either charge. Because the Compton contribution to the 
ep amplitude is linear and the Bethe-Heitler contribution quadratic in the lepton 
charge, the interference term is projected out in the difference da{e~^p) —da{e~p) 
of cross sections, whereas it is absent in their sum. Both the "beam-spin" asym- 
metry 

do+{<'^p) - da^[e-p) 
d(T+(e~p) -|- da-.{e~p) ' 

where ± denote positive (negative) beam helicities, and the "beam-charge" asym- 
metry 

da{e+p) - da{e-p) 

da{e+p) + da{e-p) ^ ' 

have been observed (47-49) . Fig. 10 shows some of the results. 

Hard exclusive meson production, ^*p Mp, is another process that gives 
access to GPDs, and much activity has gone into this direction as well (42,50). 
Both DVCS and exclusive meson production have their practical advantages and 
disadvantages. Real photon production is cleaner, but the price to be paid is an 
additional power of acm- Meson production may be easier to detect; however, 
its amplitude is suppressed relatively by a power 1/Q. The importance of using 
nucleon polarization in off-forward reactions is well established. There have also 
been first studies for DVCS off nuclei (51). 

Practical problems are the fact that GPDs depend on three variables (plus a 
scale in which they evolve), and that they appear in complicated convolutions 
with the partonic hard-scattering kernels. We are still far from the quantita- 
tive experimental surveys of DVCS and related processes that would allow us to 
work backwards to new insights into off-diagonal matrix elements and angular 
momentum. Nevertheless, a direction for the field has been set. 
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Figure 10: Data for the beam charge asymmetry in DVCS from HERMES (47) 
(left) and for the beam spin asymmetry from CLAS (49) (right), as functions of 
the azimuthal angle (j). For the definitions of these asymmetries, see text. 



2.2.3 Transverse polarization 

In addition to the unpolarized and the helicity-dependent distributions, there is 
a third set of twist-2 parton distributions, namely transversity (52). In anal- 
ogy with Equation 8 these distributions measure the net number (parallel minus 
antiparallel) of partons with transverse polarization in a transversely polarized 
nucleon: 

6q{x) = qlix) - g|(x) . (19) 

In a helicity basis (52), transversity corresponds to an interference of an amplitude 
in which a helicity- + quark emerges from a helicity- + nucleon, but is returned 
as a quark of negative helicity into a nucleon of negative helicity. This helicity- 
flip structure makes transversity a probe of chiral symmetry breaking in QCD 
(53). Perturbative-QCD interactions preserve chirality, and so the helicity flip 
must primarily come from soft non-perturbative interactions for which chiral 
symmetry is broken. The required helicity flip also precludes a gluon transversity 
distribution at leading twist (52). 

Measurements of transversity are not straightforward. Again the fact that 
perturbative interactions in the Standard Model do not change chirality (or, 
for massless quarks, helicity) means that inclusive DIS is not useful. Collins, 
however, showed (54) that properties of fragmentation might be exploited to 
obtain a "transversity polarimeter" : a pion produced in fragmentation will have 
some transverse momentum with respect to the fragmenting parent quark. There 
may then be a correlation of the form iST {P-k ^k±) among the transverse spin St 
of the fragmenting quark, the pion momentum P^, and the transverse momentum 
kj_ of the quark relative to the pion. The fragmentation function associated with 
this correlation is denoted as H^'^iz), the Collins function. If non- vanishing, the 
Collins function makes a leading-power (54-56) contribution to the single-spin 
asymmetry A±_ in the reaction ep^ — > eirX: 

oc \St\ sin(</. + <l,s) J2 el5q{x)Hi^\z) , (20) 

9 
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where (p {4's) is the angle between the lepton plane and the (7*7r) plane (and the 
transverse target spin). As shown in Equation 20, this asymmetry would then 
allow access to transversity. 

If "intrinsic" transverse momentum in the fragmentation process can play a 
crucial role in the asymmetry for ep^ — > eirX, a natural question is whether k± 
in the initial state can be relevant as well. Sivers suggested (57) that the k± dis- 
tribution of a quark in a transversely polarized hadron could have an azimuthal 
asymmetry, St {P x k±). It was realized (58,59) that the Wilson lines in the 
operators defining the Sivers function, required by gauge invariance, are crucial 
for the function to be non-vanishing. This intriguing discovery has been one of 
the most important theoretical developments in QCD spin physics in the past 
years. Another important aspect of the Sivers function is that it arises as an 
interference of wave functions with angular momenta Jz = ±1/2 and hence con- 
tains information on parton orbital angular momentum (58,60), complementary 
to that obtainable from DVCS. 

Model calculations and phenomenological studies of the Sivers functions /{^^ 
have been presented (61). It makes a contribution to ep^ eirX (55), 

oc \St\ sm{<p - CPs) E 4 firi^) D^z) . (21) 

This is in competition with the Collins function contribution. Equation 20; how- 
ever, the azimuthal angular dependence is disccrnibly different. HERMES has 
completed a run with transverse polarization and performed an extraction of the 
contributions from the Sivers and Collins effects (62) . There are also first results 
from COMPASS (63). First independent information on the Collins functions is 
now coming from Belle measurements in e+e" annihilation (64). The Collins and 
Sivers functions are also likely involved (65) in explanations of experimental ob- 
servations of very large single-transverse spin asymmetries in pp scattering (66) , 
where none were expected. It was pointed out (59, 67) that comparisons of DIS 
results and results from p^p scattering at RHIC will be particularly interesting: 
from the properties of the Wilson lines it follows that the Sivers functions vio- 
late universality of the distribution functions. For example, the Sivers functions 
relevant in DIS and in the Drell-Yan process should have opposite sign. This is 
a striking prediction awaiting experimental testing. 

2.3 Nuclear modifications 

The nucleus is traditionally described as a collection of weakly bound nucleons 
confined in a potential created by their mutual interaction. It came as a surprise 
when the EMC experiment (68) uncovered a systematic nuclear dependence to 
the nuclear structure function F2'{x,Q^) in iron relative to that for Deuterium 
because the effect was as much as 20% for x ~ 0.5. This is significantly larger than 
the effect (< 5%) due to the natural scale for nuclear effects given by the ratio 
of the binding energy per nucleon to the nucleon mass. Several dedicated fixed 
target experiments (69-71) confirmed the existence of the nuclear dependence 
observed by the EMC albeit with significant modifications of the original EMC 
results at small x. The upper part of Figure 11 shows an idealized version of the 
nuclear modification of the relative structure functions per nucleon. It is 2/ A 
times the ratio of a measured nuclear structure function of nucleus A to that for 
Deuterium. The rise at the largest values of x is ascribed to the nucleons' Fermi 
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momentum. The region above x > 0.2 is referred to as the EMC effect region. 
When X < 0.05, the nuclear ratio drops below one and the region is referred to 
as the nuclear shadowing region, whereas the region with the slight enhancement 
in between the shadowing and EMC effect regions is called the anti-shadowing 
region. The lower part of Figure 11 presents a sample of high precision data of 
ratios of structure functions over a broad range in A, x and Q^. We shall now 
discuss what is known about these regions, focusing in particular on the EMC 
effect and nuclear shadowing regions. 




Figure 11: Upper: An idealized depiction of the ratio of the structure function of 
a nucleus, F2^(x,Q^) per nucleon to F2(x,Q^) of deuterium. Lower: Measured 
F2{x, Q'^) structure functions for C, Ca, and Xe relative to Deuterium. From (72). 



2.3.1 The EMC effect 

A review of the DIS data and various interpretations of the EMC effect since its 
discovery in the early 1980's can be found in Ref. (73). A common interpreta- 
tion of the EMC effect is based on models where inter-nucleon interactions at 
a wide range of inter-nucleon distances are mediated by meson exchanges. The 
traditional theory (74, 75) of nuclear interactions predicts a net increase in the 
distribution of virtual pions with increasing nuclear density relative to that of free 
nucleons. This is because meson interactions are attractive in nuclei. In these 
models, nuclear pions may carry about 5% of the total momentum to fit the EMC 
effect at X ~ 0.3. Each pion carries a light-cone fraction of about 0.2-0.3 of that 
for a nucleon. Sea anti-quarks belonging to these nuclear pions may scatter off 
a hard probe. Hence the predicted enhancement of the nuclear sea of 10% to 
15% for X ~ 0.1-0.2 and for A > 40. The conventional view of nuclear binding is 
challenged by the constancy with A of the anti-quark distribution extracted from 
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the production of Drell-Yan pairs in proton-nucleus collisions at Fermilab (76). 
These data are shown in Fig. 12. No enhancement was observed at the level of 1% 
accuracy in the Drell-Yan experiments. The Drell-Yan data was also compared 
with and showed good agreement with the DIS EMC data for the F2 ratio of Tin 
to Deuterium. 
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Figure 12: The ratio of the anti-quark distribution per nucleon in several nuclei 
relative to Deuterium. Data are shown from a Drell-Yan experiment (76) and 
compared to theoretical predictions. Also shown is the ratio for Tungsten to 
Deuterium compared to DIS data from the EMC experiment (68) for the F2 
ratio of Tin to Deuterium. 



Furthermore, first results (77) from the Jefferson Laboratory (TJNAF) exper- 
iment E91-003 indicate that there is no significant pion excess in the A{e, e', B) 
reaction. (It has however been pointed out (78) that parameters of pion interac- 
tions in nuclei can be readily adjusted to reduce the pion excess to conform with 
the Drell-Yan data.) In addition, the energy excitation for the residual nuclear 
system also reduces the contribution of pions to the nuclear parton densities (79). 
Thus all of these observations suggest that pions may not contribute significantly 
to in the EMC region. 

The chiral quark-soliton model (33) is a phenomenologically successful model 
that for instance explains the difference in the anti-quark up and down distribu- 
tions as a function of x (80). Interestingly, it has been shown recently (80) to 
simultaneously provide a good description of both the EMC effect and the ratio of 
anti-quark distributions from Drell-Yan pairs. Recently, it has been argued that a 
key feature of the EMC effect, the factorization of the x and A dependence of the 
EMC ratio, can be understood in a model independent way in an effective field 
theory approach (81). Several joint leading order QCD analyses of the nuclear 
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DIS and Drell-Yan data combined with the apphcation of the baryon charge and 
momentum sum rules (82—85) provide further information on the nuclear effects 
on parton densities in this kinematic region. These analyses indicate that the 
valence quark distribution in nuclei is enhanced at x ~ 0.1-0.2. Gluons in nuclei 
carry practically the same fraction of the momentum (within 1%) as in a free 
nucleon. If one assumes that gluon shadowing is similar to that for quarks, these 
analyses predict a significant enhancement of the gluon distribution in nuclei at 
X ~ 0.1-0.2 (86). A recent next-to-leading order (NLO) analysis of nuclear par- 
ton distributions (87) however finds that this gluon "anti-shadowing" is much 
smaller than in the LO analysis. 

2.3.2 Nuclear shadowing 

Nuclear shadowing is the phenomenon, shown in Fig. 11, where the ratio of the 
nuclear electromagnetic structure function relative to A/2 times the Deuteron 
electromagnetic structure functions F2' is less than unity for x < 0.05. Shadowing 
is greater for decreasing x and with increasing nuclear size. For moderately small 
X, shadowing is observed to decrease slowly with increasing Q^. Unfortunately, 
because x and are inversely correlated for fixed energies, much of the very 
small X data (x < 10~^) is at very low values of Q-^ < 1 GeV^. In addition, as 
results (72) from the fixed target E665 experiment at Fcrmilab and the New Muon 
Collaboration (NMC) experiment at CERN shown in Fig. 11(b) suggest, good 
quality data exists only for a; > 4x 10~^. At high Q^, the shadowing of can be 
interpreted in terms of shadowing of quark and anti-quark distributions in nuclei 
at small x. Information on quark shadowing can also be obtained from proton- 
nucleus Drell-Yan experiments (88) and from neutrino-nucleus experiments-most 
recently from NuTeV at Fermilab (89). 

The phenomenon of shadowing has different interpretations depending on the 
frame in which we consider the space-time evolution of the scattering. Consider 
for instance the rest frame of a nucleus in 7-p/A scattering. The 7p cross-section 
is only O.lmb for energies in excess of 2 GeV, corresponding to a mean- free-path 
of well over 100 fm in nuclear matter. However, although the high-energy 7A 
cross-section might be expected to be proportional to A, the observed increase 
in the cross-section is smaller than A times the 7p cross-section. This is because 
the photon can fluctuate into a qq~ pair that has a cross-section typical of the 
strong interactions (~ 20mb) and is absorbed readily (with a mean free path 
of ~ 3.5 fm). If the fluctuation persists over a length greater than the inter- 
nucleon separation distance (2 fm), its absorption shadows it from encountering 
subsequent nucleons. The coherence length of the virtual photon's fluctuation 
is ^coh. ^ l/2mNX where miv is the nucleon mass. Therefore the onset of shad- 
owing is expected and observed at x ^ 0.05. In this Gribov multiple scattering 
picture (90), there is a close relation between shadowing and diffraction. The 
so-called AGK cutting rules (91) relate the first nuclear shadowing correction 
to the cross-section for diffractively producing a final state in coherent scatter- 
ing off a nucleon (integrated over all diffractive final states). Sec Fig. 13(a) for 
an illustration of this correspondence. With these relations (and higher order 
re-scattering generalizations of these) and with the HERA diffractive DIS data 
as input, the NMC nuclear shadowing data can be reproduced as shown in the 
sample computation (92,93) in Fig. 13(b). 

In the infinite momentum frame (IMF) , shadowing arises due to gluon recombi- 
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Figure 13: Left: An illustration of the AGK rules relating shadowing corrections 
in nuclei to diffractive scattering on nucleons. Right: A calculation (92) which 
uses AGK to fit nuclear F2 data using HERA diffractive data. The two curves 
correspond to two different unitarization prescriptions. 

nation and screening in the target. When the density of partons in the transverse 
plane of the nucleus becomes very large, many body recombination and screen- 
ing effects compete against the growth in the cross-section, leading eventually to 
a saturation of the gluon density (94). In the IMF picture, one can again use 
the AGK rules we discussed previously to relate shadowing and diffraction (95), 
and the result is amenable to a partonic interpretation. The saturation regime is 
characterized by a scale Qs{x,A), called the saturation scale, which grows with 
decreasing x and increasing A. This saturation scale arises naturally in the Color 
Glass Condensate (CGC) framework which is discussed in section 3. 

A natural consequence of saturation physics is the phenomenon of geometrical 
scaling. (See for instance the discussion on geometrical scaling of HERA data in 
section 2.1.) It has been argued that the NMC DIS data also display geometrical 
scaling (96)-the evidence here albeit interesting is not compelling owing to the 
paucity of nuclear data over a wide range of x and Q^. It is widely believed that 
shadowing is a leading twist effect (97,98), but some of the IMF discussion in 
the CGC saturation framework suggests higher twist effects are important for 
because of the large gluon density (99). Constraints from non-linear 
corrections to the DGLAP framework have also been discussed recently (100). 
The available data on the dependence of shadowing are inconclusive at small 

X. 

Our empirical definition of shadowing in DIS refers to quark shadowing. Like- 
wise for quarks and anti-quarks in the Drell-Yan process in hadronic collisions. 
In DIS gluon distributions are inferred only indirectly because the virtual photon 
couples to quarks. The most precise extractions of gluon distributions thus far 
are from scaling violations of Fg^. To do this properly, one needs a wide window 
in X and Q^. In contrast to the highly precise data on nucleon gluon distribu- 
tions from HERA, our knowledge of nuclear gluon structure functions {gA{x, Q^)) 
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Figure 14: The ratio, r{x) of the gluon distributions in Sn relative to C and 
the ratio, fi{x) of their ^2(3;) structure functions (101). The box represents the 
extraction of r{x) from J/ip electro-production in the process fi+A — ^ fj,+J/ip+X. 



is nearly non-existent. This is especially so relative to our knowledge of quark 
distributions in nuclei. The most precise data on the modification of gluon distri- 
butions in nuclei come from two NMC high precision measurements of the ratio 
of the scaling violations of the structure functions of Tin (Sn) and Carbon (C). 
The experiments measure ratios /i = i^2'V-^2^ h = Qi^gi fi- The ratio 
f = gsn/gc can be determined (101) from /i and /2 and the scaling violations of 
^Deuterium (^-^jj^j^ minimal assumptions). The result for r is shown in Fig. 14. At 
small X, gluon shadowing is observed. The trend suggests that gluon shadowing 
at small x is greater than that of F2, even though the error bars are too large for 
a conclusive statement. 

At larger x of 0.1 < x < 0.2, one observes anti-shadowing of the gluon dis- 
tributions. This result from scaling violations can be compared to the ratio of 
gluon distributions extracted from inclusive J/ip production in DIS. The latter 
assumes the gluon fusion model of J/f/j-production. The results for r from the 
latter method are consistent with those from scaling violations. The large experi- 
mental uncertainties however leave the extent of anti-shadowing in doubt. Other 
measurements of scaling violations for the ratio of Ff°/F2-' showed an increase 
of the ratio with the increase of consistent with predictions (82,83). 

The limited data we have may be interpreted to suggest a provocative picture of 
nuclear parton densities in the x ~ 0.1-0.2 region, which corresponds to distances 
of ~ 1-1.5 fm, where mediiim range and short range inter-nucleon forces are ex- 
pected to be important. In this region, if the gluon and valence quark fields are 
enhanced while the sea is somewhat suppressed, as some analyses suggest, gluon- 
induced interactions between nucleons, as well as valence quark interchanges be- 
tween nucleons may contribute significantly to nuclear binding (97,101). Nuclear 
gluon distributions can also be further constrained by inclusive hadron distribu- 
tions recently measured by the RHIC experiments (102-105) in Deuteron-Gold 
scattering at s/s = 200 GeV/nucleon. These RHIC results will be discussed in 
section 3. 
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2.4 Space-time correlations in QCD 

The space-time picture of DIS processes strongly depends on the value of Bjorken 
X. An analysis of electromagnetic current correlators in DIS reveals that one 
probes the target wave function at space-time points separated by longitudinal 
distances /coh. a-nd transverse distances ~ l/Q. At large (x > 0.2), the vir- 
tual photon transforms into a strongly interacting state very close to the active 
nucleon, typically in the middle of the nucleus. If is large enough as well, 
the produced partonic state interacts weakly with the medium. At smaller x 
{x < 0.05), the longitudinal length scale /coh. exceeds the nuclear size of the 
heaviest nuclei. At sufficiently small x {x < 0.005 for the heaviest nuclei) DIS 
processes undergo several spatially separated stages. First, the virtual photon 
transforms into a quark- gluon wave packet well before the nucleus. Time dilation 
ensures that interactions amongst partons in the wave packet are frozen over large 
distances. (These can be several hundred fermis at EIC energies.) The partons 
in the wave packet interact coherently and instantaneously with the target. At 
high energies, these interactions arc cikonal in nature and do not affect the trans- 
verse size of the wave packet. Finally, the fast components of the wave packet 
transform into a hadronic final state when well past the nucleus. This interval 
could be as large as Iv j ij? where v is the energy of the virtual photon and /x < 1 
GeV is a soft hadronic scale. 

Space-time studies thus far have been limited to semi-exclusive experiments 
that investigate the phenomenon of color transparency, and more generic inclusive 
studies of quark propagation through nuclei. Both these studies involved fixed 
targets. They are briefly summarized below. 

In pQCD, color singlet objects interact weakly with a single nucleon in the 
target. Additional interactions are suppressed by inverse powers of . This 
phenomenon is called "color transparency" because the nucleus appears trans- 
parent to the color singlet projectile (106-108). At very high energies, even the 
interaction of small color singlet projectiles with nuclei can be large. In this kine- 
matic region, the phenomenon is termed "color opacity" (109,110). The earliest 
study of color transparency in DIS was a study of coherent J/ V' photo-production 
off nuclei (111). The amplitude of the process at small t (momentum transfer 
squared) is approximately proportional to the nuclear atomic number A. This 
indicates that the pair that passes through the nucleus is weakly absorbed. For 
hadronic projectiles, a similar and approximately linear A-dependence of the am- 
plitude was observed recently for coherent diff'raction of 500 GeV pions into two 
jets (112), consistent with predictions (109). 

A number of papers (113-115) predict that the onset of color transparency at 
sufficiently large will give for the coherent diffractive production of vector 
mesons 

For incoherent diffraction at sufficiently large t (> 0.1 GeV^), they predict, 

The first measurements of incoherent diffractive production of vector mesons were 
performed by the E665 collaboration at Fermilab (116). A significant increase 
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of the nuclear transparency, as reflected in the ratio R(Q^), was observed. The 
hmited luminosity and center-of-mass energy however do not provide a statisti- 
cally convincing demonstration of color transparency. In addition, the results are 
complicated by large systematic effects. 

Measurements of the inclusive hadron distribution for different final states as 
a function of the virtual photon energy u, its transverse momentum squared Q^, 
the fraction Zh of the photon energy carried by the hadron, and the nuclear size 
A, provide insight into the propagation of quarks and gluons in nuclear media. In 
addition to the time and length scales discussed previously, the "formation time" 
Th of a hadron is an additional time scale. It is in principle significantly larger 
than the production time 1/Q of a color singlet parton. If the formation time is 
large, the "pre-hadron" can multiple scatter in the nucleus, thereby broadening 
its momentum distribution, and also suffer radiative energy loss before hadroniza- 
tion. The QCD prediction for transverse momentum broadening resulting from 
multiple scattering is (for quarks) given by the expression (117) 

(Api) = ^^^xgix, Q^)pL ^ 0.5 as GeV^ . (24) 

Here, Cp = 4/3 is the color Casimir of the quark, p is the nuclear matter density 
and L is the length of matter traversed. The Drell-Yan data in Ref. (118) agree 
with this expression and with the predicted small size of the effect (empirically, 
(Apj^) ~ 0.12 GeV^ for heavy nuclei). One also observes a large difference in 
the A dependence of the transverse momentum of Drell-Yan di-muons relative 
to those from J/i/; and T production and decay (119). In the former process 
only the incident quark undergoes strong interactions, whereas in the latter, the 
produced vector mesons interact strongly as well. However, the size of the effect 
and the comparable broadening of the J/ and T (albeit the latter is appreciably 
smaller than the former) need to be better understood. The p± imbalance of di- 
jets in nuclear photo-production suggests a significantly larger pj_ broadening 
effect than in J ftp production (120). This suggests non- universal behavior of 
p± broadening effects but may also occur from a contamination of the jets by 
soft fragments. Parton p± broadening due to multiple scattering may also be 
responsible for the anomalous behavior of inclusive hadron production in hadron- 
nucleus scattering at moderate p± of a few GeV. In this case, the ratio RpA of 
inclusive hadron production in hadron-nucleus scattering to the same process on 
a nucleon, is suppressed at low p± but exceeds unity between 1 — 2 GeV. This 
"Cronin effect" (121) was discovered in proton-nucleus scattering experiments 
in the late 70's. The flavor dependence of the Cronin effect provided an early 
hint that scattering of projectile partons off gluons dominates over scattering off 
quarks (122). The Cronin effect will be discussed further in section 3 in light of 
the recent RHIC experiments on Deuteron-Gold scattering (102-105). 

The energy loss of partons due to scattering in nuclear matter is complicated 
by vacuum induced energy loss in addition to the energy loss due to scattering. 
One computation suggests that vacuum energy loss is the dominant effect (123). 
For a quark jet, the medium induced energy loss increases quadratically with 
the length, L, and is independent of the energy for E oo. For L = 5 fm, the 
asymptotic energy loss, AE , is estimated to be less than 1 GeV in a cold nuclear 
medium (117). This makes it difficult to empirically conflrm this remarkable in- 
dependence of the energy loss. DIS data are qualitatively consistent with small 
energy loss (124-126). The data indicate that the multiplicity of the leading 
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hadrons is moderately reduced (by 10%) for virtual photon energies of the order 
of 10-20 GeV for scattering off Nitrogen-14 nuclei. At higher energies, the leading 
multiplicities gradually become A-dependent, indicating absorption of the leading 
partons (125-129). 




Figure 15: Data from the HERMES experiment (124) showing the ratio of the 
inclusive hadron cross section in a nucleus relative to that in a nuclcon, plotted 
as a function a) of Zh, which is the fraction of the quark's momentum carried by 
a hadron, and b) as a function of v, the photon energy, for two different nuclei. 
Curves denote results of a "pre-hadron" scattering model (130), with differing 
formation times. 



A pQCD description of partonic energy loss in terms of modified fragmentation 
functions is claimed to describe HERMES data (131). However, at HERMES 
energies, and perhaps even at EMC energies, descriptions in terms of hadronic 
re-scattering and absorption are at least as successful (130,132). As previously 
discussed, however, the latter descriptions usually require that the color singlet 
"pre-hadrons" have a formation time Th ~ 0.5 fm (123,130,133). Fig. 15 shows 
the results from one such model as a function of Zh (the fraction of the parton 
momentum carried by a hadron) and u (the virtual photon energy) compared to 
the HERMES data. At EIC energies, a pQCD approach in terms of modified 
fragmentation functions (134, 135) should be more applicable. The results from 
these analyses will provide an important test of jet quenching in hot matter 
descriptions of the RHIC data. 

3 Scientific Opportunities with an Electron-Ion Collider 

This section will discuss the exciting scientific opportunities that will be made 

possible by the novel features of an electron-ion collider: the high luminosity, the 
possibility to do scans over a wide range in energy, polarization of the electron 
and hadronic beams, a range of light and heavy nuclear beams and, not least, the 
collider geometry of the scattering. Scientific firsts for the electron-ion collider 
will include a) the first high energy polarized electron-polarized proton collider. 
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and b) the first high energy electron-nucleus collider. 
3.1 Unpolarized ep collisions at EIC 

Unpolarized e-p collisions have been studied extensively most recently at the 
HERA collider at DESY. In the eRHIC option for an EIC, the center of mass 
energy in an e-p collision is anticipated to be y/s = 100 GeV compared to y/s 
of over 300 GeV at HERA. Although the x-Q^ reach of an EIC may not be as 
large as that of HERA, it has significant other advantages which we will itemize 
below. 

• The current design luminosity is approximately 25 times the design luminos- 
ity of HERA. Inclusive observables will be measured with great precision. 
The additional luminosity will be particularly advantageous for studying 
semi-inclusive and exclusive final states. 

• The EIC (particularly in the eRHIC version) will be able to vary the energies 

of both the electron and nuclcon beams. This will enable a first measurement 
of Fl in the small x regime. The F^ measurement is very important in 
testing QCD fits of structure functions. 

• Electron-Deuteron collisions, with tagging of spectator nucleons, will allow 
high precision studies of the flavor dependence of parton distributions. 

• An eRHIC detector proposed by Caldwell ct al. (136) would have a rapid- 
ity coverage nearly twice that of the ZEUS and HI detectors at HERA. 
This would allow the reconstruction of the event structure of hard forward 
jets with and without rapidity gaps in the final state. With this detector, 
exclusive vector meson and DVCS measurements can be performed for a 
wider range of the photon-proton center of mass energy squared VF^. It 
also permits measurements up to high where t denotes the square of the 
difference in four-momenta of the incoming and outgoing proton. These will 
enable a precise mapping of the energy dependence of final states, as well 
as open a window into the spatial distribution of partons down to very low 
impact parameters. 

We will briefly discuss the physics measurements that can either be done or 
improved upon with the above enumerated capabilities of EIC/eRHIC. For in- 
clusive measurements, Fl is clearly a first, "gold plated" measurement. Current 
QCD fits predict that Fl is very small (and in some analyses negative) at small 
X and small < 2 GeV^. An independent measurement can settle whether 
this reflects poor extrapolations of data (implying leading twist interpretations 
of data are still adequate in this regime), or whether higher twist effects are 
dominant. It will also constrain extractions of the gluon distribution because Fl 
is very sensitive to it. Another novel measurement would be that of structure 
functions in the region of large x ~ 1. These measurements can be done with 
1 fb^^ of data for up to re = 0.9 and for < 250 GeV^. This kinematic win- 
dow is completely unexplored to date. These studies can test perturbative QCD 
predictions of the helicity distribution of the valence partons in a proton (138) 
as well as the detailed pattern of SU(6) symmetry breaking (139). Moments of 
structure functions can be compared to lattice data. These should help quantify 
the influence of higher twist effects. Finally ideas such as Bloom-Gilman duality 
can be further tested in this kinematic region (140). 
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At small X, very little is understood about the quark sea. For instance, the 
origins of the u — d asymmetry and the suppression of the strange sea are not 
clear. High precision measurements of tt^, K^, Kg and open charm will help 
separate valence and sea contributions in the small-x region. We have already 
discussed Generalized Parton Distributions and DVCS measurements. The high 
luminosity, wide coverage and measurements at high \t\ will quantify efforts to 
extract a 3-D snapshot of the distribution of partons in the proton. 

3.2 Polarized ep collisions at EIC 

We expect the EIC to dramatically extend our understanding of the spin structure 
of the proton through measurements of the spin structure function gi over a 
wide range in x and Q^, of its parity- violating counterparts and 55, of gluon 
polarization AG, as well as through spin-dependent semi-inclusive measurements, 
the study of exclusive reactions, and of polarized photo-production. 

3.2.1 Inclusive spin-dependent structure functions 

We have emphasized in Sec. 2.2.1 the need for further measurements of Q"^) 
at lower x and higher Q^. A particularly important reason is that one would like 
to reduce the uncertainty in the integral Ti[Q'^) and hence in AS. However, the 
behavior of 51 (x, Q^) at small x is by itself of great interest in QCD. 

At very high energies, Rcggc theory gives guidance to the expected behavior 
of gi{x). The prediction (141) is that gi{x) is flat or even slightly vanishing at 
small X, gi{x) oc x^" with —0.5 < a < 0. It is an open question how far one 
can increase or decrease energy and still trust Regge theory. A behavior of 
the form gi ~ x~° with a < is unstable under DGLAP evolution (142, 143) 
in the sense that evolution itself will then govern the small-x behavior at higher 
Q^. Under the assumption that Regge theory expectations are realistic at some 
(low) scale Qq one then obtains "pcrtiirbative predictions" for gi{x,Q'^) at Q ^ 
Qq. The fixed-target polarized DIS data indicate that although the non-singlet 
combination g^ — is quite singular at small x, the singlet piece does appear 
to be rather flat (142, 144), so that the above reasoning applies here. It turns 
out that the leading eigenvector of small-x evolution is such that the polarized 
quark singlet distribution and gluon density become of opposite sign. For a 
sizeable positive gluon polarization, this leads to the striking feature that the 
singlet part of gi{x, Q^) is negative at small x and large (142), driven by Ag. 
Figure 16 shows this dramatic behavior for different values of = 2, 10, 20, 100 
GeV^ (145). The projected statistical uncertainties at eRHIC, corresponding to 
400 pb~^ integrated luminosity with an almost Air acceptance detector, are also 
shown. Note that 400 pb~^ can probably be collected within about one week of 
e-p running of eRHIC. Thus, at eRHIC one will be well positioned to explore the 
evolution of the spin structure function gi{x,Q'^) at small x. We note that at 
small X and toward — 0, one could also study the transition region between 
the Regge and pQCD regimes (144). 

Predictions for the small-x behavior of gi have also been obtained from a 
perturbative resummation of double-logarithms ln^'^(l/,T) appearing in the 
splitting functions (147-150) at small x in perturbative QCD. Some of these 
calculations indicate a very singular asymptotic behavior of gi{x). It has been 
shown (148), however, that subleading terms may still be very important even 
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Figure 16: Possible eRHIC data (statistical accuracy) with 250 x 10 GeV collisions 
is shown for 400 pb~^. Also shown is the evolution of gi{x,Q'^) at low x for 
different values of for a positive gluon polarization (31,146). 



far below x = 10~^. 

The neutron gi structure function could be measured at eRHIC by colliding the 
electrons with polarized Deuterons or with Helium. If additionally the hadronic 
proton fragments are tagged, a very clean and direct measurement could be per- 
formed. As can be seen from Fig. 7, information on at small x is scarce. 
The small- rc behavior of the isotriplet — is particularly interesting for the 
Bjorken sum rule and because of the steep behavior seen in the fixed-target data 
(142, 144). It is estimated (151) that an accuracy of the order of 1% could be 
achieved for the Bjorken sum rule in a running time of about one month. One 
can also turn this argument around and use the accurate measurement of the 
non-singlet spin structure function and its evolution to determine the value of 
the strong coupling constant as(Q^)(31, 146). This has been tried, and the value 
one gets from this exercise is comparable to the world average for the strong 
coupling constant. It is expected that if precision low-x data from the EIC is 
available and the above mentioned non-singlet structure functions are measured 
along with their evolutions, this may result in the most accurate value of the 
strong coupling constant as(Q^). 

Because of eRHIC's high energy, very large can be reached. Here, the DIS 
process proceeds not only via photon exchange; also the W and Z contribute sig- 
nificantly. Equation 2 shows that in this case new structure functions arising from 
parity violation contribute to the DIS cross section. These structure functions 
contain very rich additional information on parton distributions (20, 152, 153). 
As an example, let us consider charged-current (CC) interactions. Events in the 
case of W exchange are characterized by a large transverse momentum imbal- 
ance caused by the inability to detect neutrinos from the event. The charge of 
the W boson is dictated by that of the lepton beam used in the collision. For 
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W exchange one then has for the structure functions gi and in Equation 6: 

{x) = Au{x) + Ad{x) + As{x) , gf~ {x) = Au{x) - Ad{x) - As{x) . (25) 

These appear in the double-spin asymmetry as defined in Ref. (20), where the 
asymmetry can be expressed in terms of structure functions as 

W- __ 2bg\^ ' + agl^ ' 

Here a = 2{y^ — 2y + 2), b = y{2 — y) and F3 is the unpolarized parity- violating 
structure function of Equation 2. Note that the typical scale in the parton densi- 
ties is My/ here. Availability of polarized neutrons and positrons is particularly 
desirable. For example, one finds at lowest order: 

gW-,v_w^,P = An, -Ad, (27) 

gr'^ + aY'^" = Au, + Ad, (28) 
gW+,P_gW-,n ^ -[Au + Au-Ad-Ad] . (29) 

The last of these relations gives, after integration over all x and taking into 
account the first-order QCD correction (152), 
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/ dx 
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^~^)9A, (30) 



equally fundamental as the Bjorken sum rule. 

A Monte Carlo study, including the detector effects, has shown that the mea- 
surement of the asymmetry in Equation 26 and the parity violating spin struc- 
ture functions is feasible at eRHIC. Figure 17 shows simulations (154) for the 
asymmetry and the structure function g^ for CC events with an electron beam. 
The luminosity was assumed to be 2 fb~^. The simulated data shown are for 
> 225 GeV^. Similar estimates exist for Measuring this asymmetry 

would require a positron beam. The curves in the figure use the polarized parton 
distributions of (155). It was assumed that the unpolarized structure functions 
will have been measured well by HERA by the time this measurement would be 
performed at eRHIC. Standard assumptions used by the HI collaboration about 
the scattered electrons for good detection were applied. The results shown could 
be obtained (taking into account machine and detector inefficiencies) in a little 
over one month with the eRHIC luminosity. It is possible that only one or both of 
the electron-proton and positron-proton collisions could be performed, depending 
on which design of the accelerator is finally chosen (see Section 4). 



3.2.2 Semi-inclusive measurements 

As we discussed in Sec. 2.2.1, significant insights into the nucleon's spin and fla- 
vor structure can be gained from semi-inclusive scattering ep — ^ ehX. Knowledge 
of the identity of the produced hadrons h allows separation of the contributions 
from the different quark flavors. In fixed target experiments, the so-called current 
hadrons are at forward angles in the laboratory frame. This region is difficult to 
instrument adequately, especially if the luminosity is increased to gain significant 
statistical accuracy. A polarized ep collider has the ideal geometry to overcome 
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Figure 17: Simulations (154) for the spin asymmetry of Equation 26 and 

the structure function as functions of logiQ(x). 



these shortfahs. The colUder kinematics open up the final state into a large solid 
angle in the laboratory which, using an appropriately designed detector, allows 
complete identification of the hadronic final state both in the current and target 
kinematic regions of fragmentation phase space. At eRHIC energies the current 
and target kinematics are well separated and may be individually studied. At 
eRHIC higher will be available than in the fixed-target experiments, making 
the observed spin asymmetries less prone to higher- twist effects, and the inter- 
pretation cleaner. 

Figure 18 shows simulations (156) of the precision with which one could mea- 
sure the polarized quark and antiquark distributions at the EIC. The events were 
produced using the DIS generator LEPTO. The plotted uncertainties are statis- 
tical only. The simulation was based on an integrated luminosity of 1 fb^^ for 5 
GeV electrons on 50 GeV protons, with both beams polarized to 70%. Inclusive 
and semi-inclusive asymmetries were analyzed using the leading order "purity" 
method developed by the SMC (34) and Hermes (35) collaborations. Excellent 
precision for Aq/q can be obtained down to x ~ 0.001. The measured average 
values vary as usual per x bin; they are in the range = 1.1 GeV^ at the 
lowest x to ~ 40 GeV^ at high x. With proton beams, one has greater sensi- 
tivity to up quarks than to down quarks. Excellent precision for the down quark 
polarizations could be obtained by using deuteron or helium beams. 

With identified kaons, and if the up and down quark distributions are known 
sufficiently well, one will have a very good possibility to determine the strange 
quark polarization. As we discussed in Sec. 2.2.1, As{x) is one of the most inter- 
esting quantities in nucleon spin structure. On the right-hand-side of Figure 18, 
we show results expected for As{x) as extracted from the spin asymmetries for 
production. As in the previous figure, only statistical uncertainties are in- 
dicated. The results are compared with the precision available in the Hermes 
experiment. 

There is also much interest in QCD in more refined semi-inclusive measure- 
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Figure 18: Left: projected precision of eRHIC measurements of the polarized 
quark and antiquark distributions (156). Right: expected statistical accuracy 
of As{x) from spin asymmetries for semi-inclusive measurements for 1 fb~^ 
luminosity operation of eRHIC, and comparison with the statistical accuracy of 
the corresponding Hermes measurements. 



ments. For example, the transverse momentum of the observed hadron may be 
observed. Here, interesting azimuthal-angle dependences arise at leading twist 
(56, 157), as we discussed in Subsec. 2.2.3. At small transverse momenta, re- 
summations of large Sudakov logarithms are required (158). Measurements at 
eRHIC would extend previous results from HERA (159) and be a testing ground 
for detailed studies in perturbative QCD. 

3.2.3 Measurements of the polarized gluon distribution Ag{x,Q'^) 

One may extract A.g from scaling violations of the structure function gi{x,Q'^). 
Figure 8 shows that indeed some initial information on Ag{x,Q'^) has been ob- 
tained in this way, albeit with very poor accuracy. The uncertainty of the integral 
of Ag is probably about 100% at the moment (31). Measurements at RHIC will 
vastly improve on this. eRHIC will offer independent and complementary in- 
formation. Thanks to the large lever arm in Q^, and to the low x that can be 
reached, scaling violations alone will constrain Ag{x,Q'^) and its integral much 
better. Studies (160) indicate for example that the total uncertainty on the inte- 
gral of AG could be reduced to about 5 — 10% by measurements at eRHIC with 
integrated luminosity of 12 fb~^ (~ 2 — 3 years of eRHIC operation). 

Lepton-nucleon scattering also offers direct ways of accessing gluon polariza- 
tion. Here one makes use of the photon-gluon fusion (PGF) process, for which the 
gluon appears at leading order. Charm production is one particularly interesting 
channel (161-163). It was also proposed (162-164) to use jet pairs, produced in 
the reaction ^*g — > qq, for a determination of Ag. This process competes with 
the QCD Compton process, 7*g qg. Feynman diagrams for these processes 
are shown in Figure 19. 

In the unpolarized case, dijet production has successfully been used at HERA 
to constrain the gluon density (165). Dedicated studies have been performed for 
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Photon Gluon Fusion QCD — Compton 



Figure 19: Feynman diagrams for the photon-gluon fusion and the QCD Compton 
processes. 
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Figure 20: The statistical precision of xAg from di-jets in LO for eRHIC, for two 
different luminosities, with predictions for sets A and C of the polarized parton 
densities of Ref. (155). 



dijet production in polarized collisions at eRHIC (166), using the MEPJET (167) 
generator. The two jets were required to have transverse momenta > 3 GeV, 
pseudorapidities —3.5 < < 4, and invariant mass sjj > 100 GeV^. A An 
detector coverage was assumed. The results for the reconstructed Ag{x) are 
shown in Figure 20, assuming luminosities of 1 fb~^ (left) and 200 pb~^ (right). 
The best probe would be in the region 0.02 < x < 0.1; at higher x, the QCD 
Compton process becomes dominant. This region is indicated by the shaded 
areas in the figure. The region 0.02 < x < 0.1 is similar to that probed at RHIC. 
Measurements at eRHIC would thus allow an independent determination of A^f 
in a complementary physics environment. 

Eventually data for the scaling violations in gi{x, Q^) and for dijet production 
in DIS will be analyzed jointly. Such a combined analysis would determine the 
gluon distribution with yet smaller uncertainties. A first preliminary study for 
eRHIC (168), following the lines of (169), indeed confirms this. 

3.2.4 Exploring the partonic structure of polarized photons 

In the photoproduction limit, when the virtuality of the intermediate photon is 
small, the ep cross-section can be approximated by a product of a photon flux and 
an interaction cross section of the real photon with the proton. Measurements 
at HERA in the photoproduction limit have led to a significant improvement in 
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our knowledge of the hadronic structure of the photon. 

The structure of the photon manifests itself in so-called "resolved" contribu- 
tions to cross sections. We show this in Figure 21 for the case of photoproduction 
of hadrons. On the left, the photon participates itself in the hard scattering, 
through "direct" contributions. On the right, the photon behaves like a hadron. 
This possibility occurs because of (perturbative) short-time fluctuations of the 
photon into qq pairs and gluons, and because of (non-perturbative) fluctuations 
into vector mesons p.,4>->^ with the same quantum numbers (170). The resolved 
contributions have been firmly established by experiments in e^e" annihilation 
and ep scattering (171). 




(a) (b) 

Figure 21: Generic direct (a) and resolved (b) photon contributions to the process 
Ip I'HX. 



A unique application of eRHIC would be to study the parton distributions of 
polarized quasi-real photons, defined as (172, 173) 

/:^P{x) = fl+{x)-fmx), (31) 

where f]^ (/-^) denotes the density of a parton / = u,d, s, . . . , g with positive 
(negative) helicity in a photon with positive helicity. The A/'''(x) give information 
on the spin structure of the photon; they are completely unmeasured so far. 

Figure 22 shows samples from studies (174, 175) for observables at eRHIC that 
would give information on the Af^{x). Two models for the A/'^(x) were consid- 
ered (173), one with a strong polarization of partons in the photon ("maximal" 
set), the other with practically unpolarizcd partons ("minimal" set). On the 
left, we show the double-spin asymmetry for photoproduction of high-pT pions, 
as a function of the pion's pseudorapidity r/iab in the eRHIC laboratory frame. 
The advantage of this observable is that for negative ?7iab, in the proton back- 
ward region, the photon mostly interacts "directly", via the process 75 — > qq, 
whereas its partonic content becomes visible at positive ryiab- This may be seen 
from the figure, for which we have also used two different sets of polarized parton 
distributions of the proton (26), mainly differing in Ag{x). 

The right part of Figure 22 shows predictions for the spin asymmetry in dijet 
photoproduction at eRHIC. If one assumes the jets to be produced by a 2 2 
partonic hard scattering, the jet observables determine the momentum fractions 
Xp^-y of the partons in the proton and the photon. Selecting events with x-y < 1, 
one therefore directly extracts the "resolved" -photon contribution. At higher 
orders, this picture is somewhat diluted, but remains qualitatively intact. Such 
measurements of dijet photoproduction cross sections at HERA (176) have been 
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particularly successful in providing information on photon structure. This makes 
the spin asymmetry a good candidate for learning about the Af^ at eRHIC. In 
the figure we show results for the asymmetry in three different bins of Xy. One 
can see that with 1 fb~^ luminosity one should be able at eRHIC to establish 
the existence of polarized resolved-photon contributions, and distinguish between 
our "maximal" and "minimal" photon scenarios. For a first exploration one could 
also use the approach of "effective" parton densities considered in (174,176-178). 
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Figure 22: Left: spin asymmetry for vr photoproduction in NLO QCD for two 
sets of polarized photon densities and two different choices of spin-dependent 
proton distributions. The error bars indicate the statistical accuracy anticipated 
for eRHIC assuming an integrated luminosity of lfb~^. Right: spin asymmetry 
for dijet production as a function of the jet transverse momentum, in three bins 
of the photon momentum fraction x^. 



We finally note that measurements of the polarized total photoproduction cross 
section at high energies would also give new valuable information on the high- 
energy contribution to the Drell-Hearn-Gerasimov sum rule (144). The latter 
relates the total cross sections with photon-proton angular momentum 3/2 and 
1/2 to the anomalous magnetic moments of the nucleon (179): 

°° dv I . . . , 1 2Tr^a ^ f 204.5 /xb p 

o^iiiy) - ^i/2(^) = = I 232.8 /.b n ' (^2) 



where on the right we have given the numerical values of the sum rule. Currently, 
the experimental result for the proton is a few percent high, and the one for 
the neutron about 20% low (180). There is practically no information on the 
contribution to the sum rule from photon energies > 3 GeV; estimates based 
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on Regge theory indicate that it is possible that a substantial part comes from this 
region. Measurements at eRHIC could give definitive answers here. The HI and 
ZEUS detectors at DESY routinely take data using electron taggers situated in 
the beam pipe 6-44 meters away from the end of the detectors. They detect the 
scattered electrons from events at very low and scattering angles. If electron 
taggers were included in eRHIC, similar measurements could be performed. The 
range of such measurements at eRHIC is estimated to be 10~^ — 10~^ GeV^. 

3.2.5 Hard exclusive processes 

As we have discussed in Subsection 2.2.2, generalized parton distributions (GPDs) 
are fundamental elements of nucleon structure. They contain both the parton 
distributions and the nucleon form factors as limiting cases, and they provide 
information on the spatial distribution of partons in the transverse plane. GPDs 
allow the description of exclusive processes at large Q^, among them DVCS. It 
is hoped that eventually these reactions will provide information on the total 
angular momenta carried by partons in the proton. Sec for example Eq. (15). 

The experimental requirements for a complete investigation of GPDs are formi- 
dable. Many different processes need to be investigated at very high luminosities, 
at large enough Q^, with polarization, and with suitable resolution to determine 
reliably the hadronic final state. The main difficulty, however, for experimental 
measurements of exclusive reactions is detecting the scattered proton. If the pro- 
ton is not detected, a "missing-mass" analysis has to be performed. In case of the 
DVCS reaction, there may be a significant contribution from the Bethe-Heitler 
process. The amplitude for the Bethe-Heitler process is known and, as we dis- 
cussed in Subsection 2.2.2, one may construct beam-spin and charge asymmetries 
to partly eliminate the Bethe-Heitler contribution. Early detector design studies 
have been performed for the EIC (181). These studies indicate that the accep- 
tance can be significantly increased by adding stations of Silicon-strip-based Ro- 
man Pot Detectors away from the central detector in a HERA-like configuration. 
The detector recently proposed for low- a; and low-Q^ studies at the EIC (136) 
(for details, see Section 5) may also be of significant use to measure the scattered 
proton. Further studies are underway and will proceed along with iterations of 
the design of the interaction region and of the beam line. 

Although more detailed studies need yet to be performed, we anticipate that 
the EIC would provide excellent possibilities for studying GPDs. Measurements 
at the collider will complement those now underway at fixed target experiments 
and planned with the 12-GeV upgrade at the Jefferson Laboratory (137). 

3.3 Exploring the nucleus with an electron-ion collider 

In this section, we discuss the scientific opportunities available with the EIC in 
DIS off nuclei. At very high energies, the correct degrees of freedom to describe 
the structure of nuclei are quarks and gluons. The current understanding of 
partonic structure is just sufficient to suggest that their behavior is non-trivial. 
The situation is reminiscent of Quantum Electrodynamics. The rich science of 
condensed matter physics took a long time to develop even though the nature of 
the interaction was well understood. Very little is known about the condensed 
matter, many-body properties of QCD, particularly at high energies. There are 
sound reasons based in QCD to believe that partons exhibit remarkable collec- 
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tive phenomena at high energies. Because the EIC wih be the first electron-ion 
collider, we will be entering a terra incognita in our understanding of the prop- 
erties of quarks and gluons in nuclei. The range in x and and the luminosity 
will be greater than at any previous fixed target DIS experiment. Further, the 
collider environment is ideal for studying semi-inclusive and exclusive processes. 
Finally, it is expected that a wide range of particle species and beam energies 
will be available to stiidy carefully the systematic variation of a wide range of 
observables with target size and energy. 

We will begin our discussion in this section by discussing inclusive "bread 
and butter" observables such as the inclusive nuclear quark and gluon structure 
functions. As wc observed previously, very little is known about nuclear structure 
functions at small x and 3> Aq^j-, ~ 0.04 GeV^. This is especially true of 
the nuclear gluon distribution. We will discuss the very significant contributions 
that the EIC can make in rectifying this situation. A first will be a reliable 
extraction of the longitudinal structure function at small x. Much progress has 
been made recently in defining universal diffractive structure functions (50, 182, 
183). These structure functions can be measured in nuclei for the first time. 
Generalized parton distributions will help provide a three-dimensional snapshot 
of the distribution of partons in the nucleus (43). 

We will discuss the properties of partons in a nuclear medium and the experi- 
mental observables that will enable us to tease out their properties. These include 
nuclear fragmentation functions that contain valuable information on hadroniza- 
tion in a nuclear environment. The momentum distributions of hadronic final 
states as functions of .x, Q^, and the fraction of the parton energy carried by a 
hadron also provide insight into dynamical effects such as parton energy loss in 
the nuclear medium. 

A consequence of small x evolution in QCD is the phenomenon of parton 
saturation (94). This arises from the competition between attractive Brcms- 
strahlung (184) and repulsive screening and recombination (many body) effects 
(95), which results in a phase space density of partons of order 1/as- At such 
high parton densities, the partons in the wavefunction form a Color Glass Con- 
densate (CGC) for reasons wc will discuss later (185). The CGC is an effective 
theory describing the remarkable universal properties of partons at high energies. 
It provides an organizing principle for thinking about high energy scattering and 
has important ramifications for colliders. The evolution of multi-parton correla- 
tions predicted by the CGC can be studied with high precision in lepton-nucleus 
collisions. 

Experimental observables measured at the EIC can be compared and con- 
trasted with observables extracted in proton/Deuteron- nucleus and nucleus- nuc- 
leus scattering experiments at RHIC and LHC. The kinematic reach of the EIC 
will significantly overlap with these experiments. Measurements of parton struc- 
ture functions and multi-parton correlations in the nuclear wave function will 
provide a deeper understanding of the initial conditions for the formation of a 
quark gluon plasma (QGP). Final state interactions in heavy ion collisions such 
as the energy loss of leading hadrons in hot matter (often termed "jet quench- 
ing") arc considered strong indicators of the formation of the QGP. The EIC will 
provide benchmark results for cold nuclear matter which will help quantify en- 
ergy loss in hot matter. Finally, recent results on inclusive hadron production in 
RHIC D-Au collisions at 200 GeV/nucleon show hints of the high parton density 
effects predicted by the CGC. We will discuss these and consider the similarities 
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and differences between a p/D-A and an e-A collider. 



3.3.1 Nuclear parton distributions 

The range of the EIC in x and was discussed previously (see Fig. 2). It is 
significantly larger than for the previous fixed target experiments. The projected 
statistical accuracy, per inverse picobarn of data, of a measurement of the ratio 
gf^-i versus x at the EIC relative to data from previous NMC measurements and 
a hypothetical future e-A collider at HERA energies is shown in Fig. 23. Here 
R denotes the ratio of nuclear structure functions, R = F2' / F2 ■ As discussed 
previously, the logarithmic derivative with of this ratio can be used to ex- 
tract the nuclear gluon distribution. The EIC is projected to have an integrated 
luminosity of several hundred pb^^ for large nuclei, so one can anticipate high 
precision measurements of nuclear structure functions at small x. In particular, 
because the energy of the colliding beams can be varied, the nuclear longitudinal 
structure function can be measured for the first time at small x. At small x and 
large Q^, it is directly proportional to the gluon distribution. At smaller values 
of Q^, it may be more sensitive to higher twist effects than F2 (11)- 
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Figure 23: The projected statistical accuracy of qIJq^ &s a function of x for 
an integrated luminosity of 1 pb~^ at the EIC (186). The simulated data are 
compared to previous data from the NMC and to data from a hypothetical e-A 
collider at HERA energies. 



Measurements of nuclear structure functions in the low x kinematic region will 
test the predictions of the QCD evolution equations in this kinematic region. The 
results of QCD evolution with depend on input from the structure functions at 
smaller values of for a range of x values. The data on these is scarce for nuclei. 
These results are therefore very sensitive to models of the small x behavior of 
structure functions at low Q^. A nice plot from Ref. (187) reproduced in Fig. 24 
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clearly illustrates the problem. Figure 24 shows results from theoretical models 
for the ratio of the gluon distribution in Lead to that in a proton as a function 
of X. Though all the models employ the same QCD evolution equations, the 
range in uncertainty is rather large at small x - about a factor of 3 at x ~ 10~^. 
Although one can try to construct better models, the definitive constraint can 
only come from experiment. 




Figure 24: Ratio of the gluon distribution in Lead to that in a proton, normal- 
ized by the number of nucleons, plotted as a function of x for a fixed = 5 
GeV2. From Ref. (187). Captions denote models-HKM (188), EKS98 (189), 
Sarcevic (190), Armesto (191), Frankfurt (93), Hijing (192). The vertical bands 
denote the accessible x regions at central rapidities at RHIC and LHC. 

The shadowing of gluon distributions shown in Fig. 24 is not understood in a 
fundamental way. We list here some relevant questions which can be addressed 
by a future electron-ion collider. 

• Is shadowing a leading twist effect, namely, is it unsuppressed by a power of 

? Most models of nuclear structure functions at small x assume this is 
the case. (For a review, see Ref. (98).) Is there a regime of x and Q^, where 
power corrections due to high parton density effects can be seen? (94,95,193) 

• What is the relation of shadowing to parton saturation? As we will discuss, 
parton saturation dynamically gives rise to a semi-hard scale in nuclei. This 
suggests that shadowing at small x can be understood in a weak coupling 
analysis. 

• Is there a minimum to the shadowing ratio for fixed and A with decreas- 
ing X? If so, is it reached faster for gluons or for quarks? 



40 



A. Deshpande et al. 



• The Gribov relation between shadowing and diffraction that we discussed 
previously is well established at low parton densities. How is it modified 
at high parton densities? The EIC can test this relation directly by mea- 
suring diffractive structure functions in ep (and e-A) and shadowing in e-A 
collisions. 

• Is shadowing universal? For instance, would gluon parton distribution func- 
tions extracted from p-A collisions at RHIC be identical to those extracted 
from e-A in the same kinematic regime? The naive assumption that this 
is the case may be false if higher twist effects are important. Later in this 
review, we will discuss the implications of the possible lack of universality 
for p-A and A- A collisions at the LHC. 

We now turn to a discussion of diffractive structure functions. At HERA, hard 
diffractive events were observed where the proton remained intact and the virtual 
photon fragmented into a hard final state producing a large rapidity gap between 
the projectile and target. A rapidity gap is a region in rapidity essentially devoid 
of particles. In pQCD, the probability of a gap is exponentially suppressed as 
a function of the gap size. At HERA though, gaps of several units in rapid- 
ity are relatively unsuppressed; one finds that roughly 10% of the cross-section 
corresponds to hard diffractive events with invariant masses Mx > 3 GeV. The 
remarkable nature of this result is transparent in the proton rest frame: a 50 TcV 
electron slams into the proton and, 10% of the time, the proton is unaffected, 
even though the interaction causes the virtual photon to fragment into a hard 
final state. 

The interesting question in diffraction is the nature of the color singlet object 
(the "Pomeron" ) within the proton that interacts with the virtual photon. This 
interaction probes, in a novel fashion, the nature of confining interactions within 
hadrons. (We will discuss later the possibility that one can study in diffractive 
events the interplay between strong fields produced by confining interactions and 
those generated by high parton densities.) In hard diffraction, because the invari- 
ant mass of the final state is large, one can reasonably ask questions about the 
quark and gluon content of the Pomeron. A diffractive structure function F^^^ 
can be defined (182,183,194), in a fashion analogous to F2, as 

d'^CTeA^eXA _ 47rag„, j _^ 

dxBjdQ^dxpdt xQ4 I y 2[l + R^''^\p,Q^,Xp,t)] 

X F^}^\p,Q\xp,t), (33) 

where, y = Q'^/sxbj, and analogously to F2, one has R^^^^ = F^^"^^ /F^^^\ 
Further, = _q2 > g, xsj = Q'^ /'iPq, Xp = q{P - P')/qP, t = {P - P'f 
and j3 = XBj/xp. Here P is the initial nuclear momentum, and P' is the net 
momentum of the fragments Y in the proton fragmentation region. Similarly, Mx 
is the invariant mass of the fragments X in the electron fragmentation region. 
An illustration of the hard diffractive event is shown in Fig. 25. 

It is more convenient in practice to measure the structure function F^^^ = 

I ^2^A^dt, where \tmin\ < \t\ < \tmax\, where \tmin\ is the minimal momen- 
tum transfer to the nucleus, and \tmax\ is the maximal momentum transfer 
to the nucleus that still ensures that the particles in the nuclear fragmenta- 
tion region Y are undetected. An interesting quantity to measure is the ra- 




Figure 25: The diagram of a process with a rapidity gap between the systems X 
and Y. The projectile nucleus is denoted here as p. Figure from Ref. (195). 



tio Ra\,A2{P,Q'^ ,Xp) = ola). ' 2' ^ • The ^l-dependence of this quantity will 

contain very useful information about the universality of the structure of the 
Pomeron. In a study for e-A collisions at HERA, it was argued that this ratio 
could be measured with high systematic and statistical accuracy (195)-the situ- 
ation for cRHlC should be at least comparable, if not better. Unlike F2 however, 
F2 is not truly universal- it cannot be applied, for instance, to predict diffractive 
cross sections in p-A scattering; it can be applied only in other lepton-nucleus 
scattering studies (50, 182). This has been confirmed by a study where diffractive 
structure functions measured at HERA were used as an input in computations 
for hard diffraction at Fermilab. The computations vastly overpredicted the Fer- 
milab data on hard diffraction (196). Some of the topics discussed here will be 
revisited in our discussion of high parton densities. 



3.3.2 Space-time evolution of partons in a nuclear environment 

The nuclear structure functions are inclusive observables and are a measure of 
the properties of the nuclear wavefunction. Less inclusive observables, which 
measure these properties in greater detail, will be discussed in the section on 
the Color Glass Condensate. In addition to studying the wave function, we are 
interested in the properties of partons as they interact with the nuclear medium. 
These are often called final state interactions to distinguish them from the initial 
state interactions in the wavefunction. Separating which effects arise from the 
wavefunction is not easy because our interpretation of initial state and some final 
state interactions may depend on the gauge in which the computations are per- 
formed (197). Isolating the two effects in experiments is difficult. A case in point 
is the study of energy loss effects on final states in p-A collisions (198). These 
effects are not easy to distinguish from shadowing effects in the wavefunction. 
Nevertheless, in the right kinematics this can be done. 

In section 2, we discussed various final state in-medium QCD processes such 
as color transparency, parton energy loss and the medium modification of frag- 
mentation functions. The EIC will enable qualitative progress in studies of the 
space-time picture of strong interactions relative to previous fixed target DIS 
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experiments. The reasons for this fohows. 

• The high luminosity of the EIC will increase by many orders of magnitude 
the current data sample of final states in DIS scattering off nuclei at high 
energies. 

• The EIC will provide a much broader range of and x, making it possi- 
ble to compare dynamics for approximately the same space-time coherence 
lengths as a function of . Fixing the coherence length of partons will allow 
one to distinguish events wherein a photon is transformed into a strongly 
interacting system either outside or inside the miclcTis. This will help isolate 
initial state interactions from those in the final state. 

• The collider geometry will enable measurements of final states currently 
impossible in fixed target kinematics. In particular, a hermetic detector 
would clearly isolate coherent processes as well as quasi-elastic processes in 
DIS off nuclei. In addition, one can study the sizes and distributions of 
rapidity gaps as a function of nuclear size and energy. These will provide a 
sensitive probe of the interplay between space-time correlations in the final 
state and in the nuclear wavefunction. 

• The detection of nucleons produced in the nuclear fragmentation region 
would make it feasible to study DIS as a function of the number of the 
nucleons involved in the interaction. In particular, it may be possible to 
study impact parameter dependence of final states, which will be important 
to understand in detail the nuclear amplification of final state effects. In 
addition, the impact parameter dependence will help distinguish geometrical 
effects from dynamical effects in event-by-event studies of final states. 

In section 2, we discussed Generalized Parton Distributions (GPDs) in the 
context of DIS scattering off nucleons. These GPDs can also be measured in DIS 
scattering off nuclei (42). The simplest system in which to study GPDs is the 
deuteron. The transition from D — > p + n in the kinematics where the neutron 
absorbs the momentum transfer in the scattering is sensitive to the GPD in the 
neutron with the proton playing the role of a spectator (199). Some preliminary 
studies have been done for heavier nuclei (51). Certain higher twist correlations 
in nuclei which scale as A^^^ are sensitive to nucleon GPDs (134). This leads 
us to a discussion of GPDs in nuclei at small x. As we will discuss, high parton 
density effects are enhanced in large nuclei. k± dependent GPDs might provide 
the right approach to study this novel regime (200). The study of nuclear GPDs 
at moderate and small a; is a very promising, albeit nascent, direction for further 
research to uncover the detailed structure of hard space-time processes in nuclear 
media. These nuclear distributions can be studied for the first time with the EIC. 

3.3.3 The Color Glass Condensate 

The Color Glass Condensate (CGC) is an effective field theory describing the 
properties of the dominant parton configurations in hadrons and nuclei at high en- 
ergies (185). The degrees of freedom are partons-which carry color charge-hence 
the "Color" in CGC. The matter behaves like a glass for the following reason. 
The kinematics of high energy scattering dictates a natural separation between 
large x and small x modes (201). The large x partons at high energies behave 
like frozen random light cone sources over time scales that are large compared to 
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the dynamical time scales associated with the small-x partons. One can therefore 
describe an effective theory where the small x partons are dynamical fields and 
the large-x partons are frozen sources (193). Under quantum evolution (202), this 
induces a stochastic coupling between the wee partons via their interaction with 
the sources. This stochastic behavior is very similar to that of a spin glass. Fi- 
nally, the Condensate in CGC arises because each of these colored configurations 
is very similar to a Bose-Einstein Condensate. The occupation number of the 
gluons can be computed to be of order 1/as, and the typical momentum of the 
partons in the configuration is peaked about a typical momentum-the saturation 
momentum Qs- These properties are further enhanced by quantum evolution in 
X. Because the occupation number is so large, by the correspondence principle 
of quantum mechanics, the small x modes can be treated as classical fields. The 
classical field retains its structure while the saturation scale, generated dynami- 
cally in the theory, grows with energy: Qs{x') > Qsix) for x' < x. The CCC is 
sometimes used interchangeably with "saturation" (203) - both refer to the same 
phenomenon, the behavior of partons at large occupation numbers. 

The Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) re- 
normalization group equations describe the properties of partons in the high den- 
sity regime (202). They form an infinite hierarchy (analogous to the Bogoliubov- 
Born-Green-Kirkwood- Young (BBGKY) hierarchy in statistical mechanics) of 
ordinary differential equations for the gluon correlators {AiA2 - ■ ■ An)Y ^ where 
Y = ln(l/a;) is the rapidity. Thus the evolution, with x, of multi-gluon (semi- 
inclusive) final states provides precise tests of these equations. The full hierarchy 
of equations arc difficult to solve^ though there have been major theoretical de- 
velopments in that direction recently (205). 

A mean field version of the JIMWLK equation, called the Balitsky-Kovchegov 
(BK) equation (206), describes the inclusive scattering of the quark-anti-quark 
dipole off the hadron in deeply inelastic scattering. In particular, the virtual 
photon-proton cross-section at small x can be written as (207, 208) 

'^T,L = j J dz\'^T,Lir±,z,Q^)faggNir_L,x), (34) 

where |V'T,Lp is the probability for a longitudinally (L) or transversely (T) po- 
larized virtual photon to split into a quark with momentum fraction z and an 
anti-quark with momentum fraction 1 — z oi the longitudinal momentum of the 
virtual photon. For the quark and anti-quark located at x± and y± respectively 
from the target, their transverse size is f± = x± — y±, and the impact parameter 
of the collision is 6 = {x± + y±)/2. The probability for this splitting is known 
exactly from QED and it is convoluted with the cross-section for the qq-pair to 
scatter off the proton. This cross-section for a dipole scattering off a target can 
be expressed as 

crqqN{x,r±) = 2 J (fb MYix,r±,b) , (35) 

where My is the imaginary part of the forward scattering amplitude. The BK 
equation (206) for this amplitude has the operator form 



Y 



dY 



as /Cbfkl <8) {My - } ■ (36) 



^For a preliminary numerical attempt, see Ref. (204). 
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Here /Cbfkl is the well known Balitsky-Fadin-Kuraev-Lipatov (BFKL) kernel (184) 
When M « I, the quadratic term is negligible and one has BFKL growth of the 
number of dipoles; when M is close to unity, the growth saturates. The approach 
to unity can be computed analytically (209). The BK equation is the simplest 
equation including both the Bremsstrahlung responsible for the rapid growth of 
amplitudes at small x as well as the repulsive many body effects that lead to a 
saturation of this growth. 

Saturation models, which incorporate key features of the CGC, explain several 
features of the HERA data. In section 2.1, we discussed the property of geomet- 
rical scaling observed at HERA which is satisfied by the LHS of Eq. 34, where it 
scaled as a function of the ratio of to the saturation scale Q^. Wc also men- 
tioned briefly a simple saturation model, the Golec-Biernat model (13), which 
captured essential features of this phenomenon in both inclusive and diffractive 
cross sections at HERA. Geometric scaling arises naturally in the Color Glass 
Condensate (210,211), and it has been studied extensively both analytically (212) 
and numerically (213-215) for the BK equation. The success of saturation mod- 
els, as discussed in section 2.1, in explaining less inclusive features of the HERA 
data is also encouraging since their essential features can be understood to fol- 
low from the BK equation. Below we will discuss the implications of mean field 
studies with the BK equation, as well as effects beyond BK. 

As mentioned previously, a very important feature of saturation is the dynam- 
ical generation of a dimcnsionful scale ^ ^QCD' which controls the running 
of the coupling at high energies: as{Ql) <C 1. From the BK equation, or more 
generally, from solutions of BFKL in the presence of an absorptive boundary 
(corresponding to a CGC-like regime of high parton densities), one can deduce 
that, for fixed coupling, has the asymptotic form = Qgexp (cY), where 
c = 4.8as and Y = ln(xo/x). Here, Qq and xq are parameters from the initial 
conditions. Pre-asymptotic Y dependent corrections can also be computed and 
are large. The behavior of Q'^ changes qualitatively when running coupling effects 
are taken into account. The state of the art is a computation of the saturation 
scale to next-to-leading order in BFKL with additional resummation of collinear 
terms that stabilize the predictions of NLO BFKL (216). One recovers the form 
Ql = (5oexp(Ay), now with small pre-asymptotic corrections, with A 0.25. 
Remarkably, this value is very close to the value extracted in the Golec-Biernat 
model from fits to the HERA data. 

Fig. 26 shows a schematic plot of the CGC and extended scaling regions in the 
x-Q^ plane. Clearly, with the wide kinematic range of the EIC, and the large 
number of available measurements-to be discussed later-one has the opportunity 
to make this plot quantitative. One can further add an additional axis for the 
atomic number to see how the kinematic reach of the CGC scales with A. In 
principle, one can also study the impact parameter dependence of the saturation 
scale in addition to the ^-dependence. 

3.3.4 Signatures of the CGC 

Inclusive signatures. Inclusive measurements include F2 and for a wide range 
of nuclei, the latter measurements being done independently for the first time. 
The data will be precise enough to extract derivatives of these with respect to 
InQ^ and In x in a wide kinematic range in x and Q^. Logarithmic derivatives of 
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Figure 26: A schematic plot of the Color Glass Condensate and extended scal- 
ing regimes in the x-Q^ plane. Here r = ln(l/x) denotes the rapidity. From 
Ref. (185). 



F2 and Fl will enable the extraction of the coefficient A of the saturation scale, 
which as discussed previously, is defined to be = Qq^^^, where Y = ln(xo/x), 
and where xq and Qq are reference values corresponding to the initial conditions 
for small-x evolution. Simulations suggest that a precise extraction of this quan- 
tity may be feasible (136). Except at asymptotic energies, A = X(Y). Predictions 
exist for "universal" pre-asymptotic y-dependent corrections to A (212). Second 
derivatives of F2 and Fl with respect to ln(xo/x) will be sensitive to these cor- 
rections. The logarithmic derivatives of F2 and Fi with Q^, especially the latter, 
will be sensitive to higher twist effects for ~ Qlix, A). The saturation scale is 
larger for smaller x and larger ^~thus deviations of predictions of CGC fits from 
DGLAP fits should systematically increase as a function of both. CGC fits have 
been shown to fit HERA data at small x (217,218). These fits can be extended 
to nuclei and compared to scaling violation data relative to DGLAP fits. The A 
dependence of the saturation scale can also be extracted from nuclear structure 
functions at small x. Again, predictions exist for the pre-asymptotic scaling of 
the saturation scale with rapidity (or x), for different A (219), that can be tested 
against the data. 

In the BK equation (mean field approximation of the CGC renormalization 
group equations), we now have a simple way to make predictions for the effects 
of high parton densities on both inclusive and diffractive (220) structure functions. 
There are now a few preliminary computations for e-A DIS in this framework (221, 
222). Much more remains to be done-in particular, comparisons with DGLAP 
for EIC kinematics and detector cuts. 

Semi-inclusive and exclusive signatures. The collider geometry of the EIC will 
greatly enhance the semi-inclusive final states in e-A relative to previous fixed tar- 
get experiments. Inclusive hadron production at p± ~ Qs should be sensitive to 
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higher twist effects for « Q'^{x,A). For the largest nuclei, these effects should 
be clearly distinguishable from DGLAP based models. Important semi-inclusive 
observables are coherent (or diffractive) and inclusive vector meson production, 
which are sensitive measures of the nuclear gluon density (110,115). Exclusive 
vector meson production was suggested by Mueller, Munier and Stasto (17) as 
a way to extract the S-matrix (and therefore the saturation scale in the Golec- 
Biernat-Wiisthoff parameterization) from the t-depcndcncc of exclusive p-mcson 
production. A similar analysis of J/ip production was performed by Guzey et 
al. (223). These studies for e-A collisions will provide an independent measure 
of the energy dependence of the saturation scale in nuclei. An extensive recent 
theoretical review of vector meson production of HERA (relevant for EIC studies 
as well) can be found in Ref. (224). 

In hard diffraction, for instance, one should be able to distinguish predictions 
based on the strong field effects of BK (or hard pomeron based approaches in 
general) from the soft pomeron physics associated with confinement. As we 
discussed previously, some saturation models predict that hard diffractive events 
will constitute 30-40% of the cross-section (225,226). These computations can 
be compared with DGLAP predictions which match soft Pomeron physics with 
hard perturbative physics. One anticipates that the latter would result in a much 
smaller fraction of the cross-section and should therefore be easily distinguishable 
from CGC based "strong field" diffraction. 

The BK renormalization group equation is not sensitive to multi-particle cor- 
relations. These are sensitive to effects such as Pomeron loops (205), although 
phenomenological consequences of these remain to be explored. These effects are 
reflected in multiplicity fluctuations and rapidity correlations over several units 
in rapidity (91,227). One anticipates quantitative studies of these will be devel- 
oped in the near future. A wide detector coverage able to resolve the detailed 
structure of events will be optimal for extracting signatures of the novel physics 
of high parton densities. 

3.3.5 Exploring the CGC in proton/Deuteron- nucleus collisions 

Although high parton density hot spots may be studied in pp collisions, they 
are notoriously hard to observe. The proton is a dilute object, except at small 
impact parameters, and one needs to tag on final states over a wide 47r cov- 
erage. Deuteron-nucleus experiments are more promising in this regard. They 
have been performed at RHIC and may be performed at LHC in the future. The 
Cronin effect discovered in the late 70's (121) predicts a hardening of the trans- 
verse momentum spectrum in proton-nucleus collisions, relative to proton-proton 
collisions at transverse momenta of order ~ 1 — 2 GeV. It disappears at much 
larger p±. A corresponding depletion is seen at low transverse momenta. The 
effect was interpreted as arising from the multiple scatterings of partons from the 
proton off partons from the nucleus (122). 

First data from RHIC on forward D-Au scattering at s/s = 200 GeV/nucleon 
demonstrate how the Cronin effect is modified with energy or, equivalently, with 
the rapidity. The x values in nuclei probed in these experiments, at ~ 2 GeV, 
range from 10~^ in the central rapidity region down to 10"*^ at very forward 
rapidities^ At central rapidities, one clearly sees a Cronin peak at ~ 1 — 2 

^It has been argued (228), however, that the forward D-A cross section in the BRAHMS 
kinematic regime receives sizable contributions also from rather large x values. 
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Figure 27: Depletion of the Cronin peak from r/ = to r/ = 3 for minimum bias 
events. From Ref. (105). 
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Figure 28: Centrality dependence of the Cronin ratio as a function of rapidity. 
From Ref. (105). 



GeV. A dramatic result obtained by the BRAHMS (105) experiment at RHIC^ is 
the rapid shrinking of the Cronin peak with rapidity shown in Fig. 27. In Fig. 28, 
the centrality dependence of the effect is shown. At central rapidities, the Cronin 
peak is enhanced in more central collisions. For forward rapidities, the trend is 
reversed: more central collisions at forward rapidities show a greater suppression 
than less central collisions! 

Parton distributions in the classical theory of the CGC exhibit the Cronin 
effect (229-231). However, unlike this classical Glauber picture (232), quantum 
evolution in the CGC shows that it breaks down completely when the X2 in the 
target is such that ln(l/x2) ~ I/els'- This is precisely the trend observed in 
the RHIC D-Au experiments (105). The rapid depletion of the Cronin effect in 
the CGC picture is due to the onset of BFKL evolution, whereas the subsequent 
saturation of this trend reflects the onset of saturation effects (233). The inversion 
of the centrality dependence can be explained as arising from the onset of BFKL 
anomalous dimensions, that is, the nuclear Bremsstrahlung spectrum changes 
from Qs/Px — ^ Qs/p±- Finally, an additional piece of evidence in support 
of the CGC picture is the broadening of azimuthal correlations (234) for which 
preliminary data now exists from the STAR collaboration (235). We note that 

"The trends seen by BRAHMS are also well corroborated by the PHOBOS, PHENIX and 
STAR experiments at RHIC in different kinematic ranges (102-104). 
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alternative explanations have been given to explain the BRAHMS data (236). 
These ideas can be tested conclusively in photon and di-lepton production in 
D-A collisions at RHIC (237) as well as by more detailed correlation studies. 

Hadronic collisions in pQCD are often interpreted within the framework of 
coUinear factorization. At high energies, k± factorization may be applicable (238) 
where the relevant quantities are "unintegrated" k± dependent parton densities. 
Strict A;_L-factorization which holds for gluon production in p-A collisions (197, 
239) is broken for quark production (240,241), for azimuthal correlations (242) 
and diffractive final states (243). For a review, see Ref. (244). These cross- 
sections can still be written in terms of A;_|_-dcpcndcnt multi-parton correlation 
functions (240) and will also appear in DIS final states (220). DIS will allow us 
to test the universality of these correlations, that is, whether such correlations 
extracted from p-A collisions can be used to compute e-A final states (230,245). 

3.3.6 The Color Glass Condensate and the Quark Gluon Plasma 

The CGC provides the initial conditions for nuclear collisions at high energies. 
The number and energy of gluons released in a heavy ion collision of identical 
nuclei can be simply expressed in terms of the saturation scale as (246-248) 

where ce ~ 0.25 and cn ~ 0.3. Here rj is the space-time rapidity. These simple 
predictions led to correct predictions for the hadron multiplicity at central ra- 
pidities in Au-Au collisions at RHIC (246,249) and for the centrality and rapidity 
dependence of hadron distributions (250). However, the failure of more detailed 
comparisons to the RHIC jet quenching data (251) and elhptic flow data (252) 
suggested that final state effects are important and significantly modify predic- 
tions based on the CGC alone. The success of hydrodynamic predictions suggests 
that matter may have thermalized to form a quark gluon plasma (253). Indeed, 
bulk features of multiplicity distributions may be described by the CGC precisely 
as a consequence of early thermalization-leading to entropy conservation (254). 
Initial-state effects will be more important in heavy ion collisions at the LHC 
because one is probing smaller x in the wave function. Measurements of satura- 
tion scales for nuclei at the EIC will independently corroborate equations such as 
Eq. 37 and therefore the picture of heavy ion collisions outlined above. Further, a 
systematic study of energy loss in cold matter will help constrain extrapolations 
of pQCD (131) used to study jet quenching in hot matter. 

3.3.7 Proton/Deuteron-nucleus versus electron-nucleus collisions as probes of 

high parton densities 

Both p/D-A and e-A collisions probe the small x region at high energies. Both 
are important to ascertain truly universal aspects of novel physics. e-A colli- 
sions, owing to the independent "lever" arm in x and Q^, as well as the simpler 
lepton-quark vertex, arc better equipped for precision measurements. For exam- 
ple, in e-A collisions, information about gluon distributions can be extracted from 
scaling violations and from photon-gluon fusion processes. In both cases, high 
precision measurements arc feasible. In p-A collisions, one can extract gluon dis- 
tributions from scaling violations in Drell-Yan and gluon-gluon and quark-gluon 
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fusion channels such as open charm and direct photon measurements respectively. 
However, for both scaling violations and fusion processes, one has more convolu- 
tions and kinematic constraints in p-A than in e-A. These limit both the precision 
and range of measurements. In Drell-Yan, in contrast to F2, clear scaling viola- 
tions in the data are very hard to see and data are limited to > 16 GeV^, 
above the J/i/: and 4>' thresholds. 

A clear difference between p/D-A and e-A collisions is in hard diffractive final 
states. At HERA, these constituted approximately 10% of the total cross section. 
At eRHIC, these may constitute 30-40% of the cross section (225,226). Also, 
factorization theorems derived for diffractive parton distributions only apply to 
lepton-hadron processes (182). Spectator interactions in p/D-A collisions will 
destroy rapidity gaps. A comparative study of p/D-A and e-A collisions thus 
has great potential for unravelling universal aspects of event structures in high 
energy QCD. 

4 Electron-Ion Collider: Accelerator Issues 

With the scientific interest in a high luminosity Icpton-ion collider gathering 
momentum during the last several years, there has been a substantial effort in 
parallel to develop a preliminary technical design for such a machine. A team of 
physicists from BNL, MIT-Bates, DESY and the Budker Institute have developed 
a realistic design (255) for a machine using RHIC, which would attain an e- 
p collision luminosity of 0.4 x 10^^ cm~^s~^ and could with minimal R&D start 
construction as soon as funding becomes available. Other more ambitious lepton- 
ion collider concepts which would use a high intensity electron linac to attain 
higher luminosity are under active consideration (255,256). This section gives an 
overview of the activities currently underway related to the accelerator design. 

The physics program described above sets clear requirements and goals for the 
lepton-ion collider to be a successful and efficient tool. These goals include: a 
sufficiently high luminosity; a significant range of beam collision energies; and po- 
larized beam (both lepton and nucleon) capability. On the other hand, to be real- 
istic, the goals should be based on the present understanding of the existing RHIC 
machine and limitations which arise from the machine itself. Realistic machine 
upgrades should be considered to overcome existing limitations and to achieve 
advanced machine parameters, hut those Tipgrades should be cost-effective. 

The intent to minimize required upgrades in the existing RHIC rings affects 
the choice of parameters and the set of goals. For example, the design assumed 
simultaneous collisions of both ion-ion and lepton-ion beams. In the main design 
line, collisions in two ion-ion interaction regions, at the "6" and "8 o'clock" 
locations, have to be allowed in parallel with electron-ion collisions. 

Taking these considerations into account, the following goals were defined for 
the accelerator design: 

• The machine should be able to provide beams in the following energy ranges: 
for the electron accelerator, 5-10 GeV polarized electrons, 10 GeV polarized 
positrons; for the ion accelerator, 50-250 GeV polarized protons, 100 GeV/u 
Gold ions. 

• Luminosity: in the 10^^ — 10^^ cm~^s~^ range for e-p collisions; in the lO^'^ — 
10^^ cm^^s"^ range for e-Au collisions 

• 70% polarization for both lepton and proton beams 
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• Longitudinal polarization in the collision point for both lepton and proton 
beams 

An additional design goal was to include the possibility of accelerating polarized 
ions, especially polarized ^He ions. 

4.1 eRHIC: ring-ring design 

The primary eRHIC design centers on a 10 GeV lepton storage ring which in- 
tersects with one of the RHIC ion beams at one of the interaction regions (IRs), 
not used by any of the ion-ion collision experiments. RHIC uses superconducting 
dipole and quadrupole magnets to maintain ion beams circulating in two rings 
on a 3834 meter circumference. The ion energy range covers 10.8 to 100 GeV/u 
for gold ions and 25 to 250 GeV for protons. There are in total 6 intersection 
points where two ion rings, Blue and Yellow, cross each other. Four of these 
intersections points are currently in use by physics experiments. 

A general layout of the ring-ring eRHIC collider is shown in Figure 29 with the 
lepton-ion collisions occurring in the "12 o'clock" interaction region. Plans have 
been made for a new detector, developed and optimized for electron-ion collision 
studies, to be constructed in that interaction region. 




Figure 29: Schematic layout of the ring-ring eRHIC collider. 



The electron beam in this design is produced by a polarized electron source 
and accelerated in a linac injector to energies of 5 to 10 GeV. To reduce the 
injector size and cost, the injector design includes recirculation arcs, so that 
the electron beam passes through the same accelerating linac sections multiple 
times. Two possible linac designs, superconducting and normal conducting, have 
been considered. The beam is accelerated by the linac to the required collision 
energy and injected into the storage ring. The electron storage ring is designed 
to be capable of electron beam storage in the energy range of 5 to 10 GeV 
with appropriate beam emittance values. It does not provide any additional 
acceleration for the beam. The electron ring should minimize depolarization 
effects in order to keep the electron beam polarization lifetime longer than the 
typical storage time of several hours. 

The injector system also includes the conversion system for positron produc- 
tion. After production the positrons are accelerated to 10 GeV energy and in- 
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jected into the storage ring similarly to the electrons. Obviously the field polari- 
ties of all ring magnets should be reversed in the positron operation mode. Unlike 
electrons, the positrons are produced unpolarized and have to be polarized using 
radiative self-polarization in the ring. Therefore, the design of the ring should 
allow for a sufficiently small self polarization time. The current ring design pro- 
vides a self polarization time of about 20 min at 10 GeV. But with polarization 
time increasing sharply as beam energy goes down the use of a polarized positron 
beam in the present design is limited to 10 GeV energy. 

The design of the eRHIC interaction region involves both accelerator and de- 
tector considerations. Figure 29 shows the electron accelerator located at the "12 
o'clock" region. Another possible location for the electron accelerator and for 
electron-ion collisions might be the "4 o'clock" region. For collisions with elec- 
trons the ion beam in the RHIC Blue ring will be used, because the Blue ring can 
operate alone, even with the other ion ring. Yellow, being down. The interaction 
region design provides for fast beam separation for electron and Blue ring ion 
bunches as well as for strong focusing at the collision point. In this design, the 
other (Yellow) ion ring makes a 3m vertical excursion around the collision region, 
avoiding collisions both with electrons and the Blue ion beam. The eRHIC inter- 
action region includes spin rotators, in both the electron and the Blue ion rings, 
to produce longitudinally polarized beams of leptons and protons at the collision 
point. 

The electron cooling system in RHIC (257, 258) is one of the essential upgrades 
required for eRHIC. The cooling is necessary to reach the luminosity goals for lep- 
ton collisions with Gold ions and low (below 150 GeV) energy protons. Electron 
cooling is considered an essential upgrade of RHIC to attain higher luminosity in 
ion-ion collisions. 

In addition, the present eRHIC design assumes a total ion beam current higher 
than that being used at present in RHIC operation. This is attained by operating 
RHIC with 360 bunches. 

The eRHIC collision luminosity is limited mainly by the maximum achievable 
beam-beam parameters and by the interaction region magnet aperture limita- 
tions. To understand this, it is most convenient to use a luminosity expression in 
terms of beam-beam parameters (4eCi) and rms angular spread in the interaction 
point (o-^,cj(,J: 

The /c = 28.15 MHz is a collision frequency, assuming 360 bunches in the ion 
ring and 120 bunches in the electron ring. The parameter K = (Jy/(Jx presents 
the ratio of beam sizes in the interaction point. One of the basic conditions which 
defines the choice of beam parameters is a requirement on equal beam sizes of 
ion and electron beams at the interaction point: Uxe = <^xi and dyf. = ay^. The 
requirement is based on the operational experience at the HERA collider and on 
the reasonable intention to minimize the amount of one beam passing through 
the strongly nonlinear field in the outside area of the counter-rotating beam. 

According to the above expression, the luminosity reaches a limiting value at 
the maximum values of beam-beam parameters, or at the beam-beam parameter 
limits. For protons (and ions) the total beam-beam parameter limit was assumed 
to be 0.02 , following the experience and observation from other proton machines 
as well as initial experience from RHIC operation. With three beam-beam inter- 
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action points, two for proton-proton and one for electron-proton collisions, the 
beam-beam parameter per interaction point should not exceed 0.007. 

For the electron (or positron) beam a limiting value of the beam-beam pa- 
rameter has been put at 0.08 for 10 GeV beam energy, following the results of 
beam-beam simulations, as well as from the experience at electron machines of 
similar energy range. Because the beam-beam limit decreases proportionally with 
the beam energy, the limiting value for 5 GeV is reduced to 0.04. 

The available magnet apertures in the interaction region also put a limit on 
the achievable luminosity. The work on the interaction region design revealed 
considerable difficulties to provide an acceptable design for collisions of round 
beams. The IR has been designed to provide low beta focusing and efficient 
separation of elliptical beams, with beam size ratio K = 1/2. The main aperture 
limitation comes from the septum magnet, which leads to the limiting values of 
cr^p = 93;urad. 

Another limitation which must be taken into account is a minimum acceptable 
value of the beta-function at the interaction point (/?*). With the proton rms 
bunch length of 20 cm, decreasing f3* well below this number results in a luminos- 
ity degradation due to the hour-glass effect. The limiting value /3* = 19 cm has 
been used for the design, which results in a luminosity reduction of only about 
12%. A bunch length of 20 cm for Au ions would be achieved with electron 
cooling. 

Tables 1 and 2 show design luminosities and beam parameters. The positron 
beam intensity is assumed to be identical to the electron beam intensity, hence the 
luminosities for collisions involving a positron beam are equivalent to electron-ion 
collision luminosities. To achieve the high luminosity in the low energy tune in 
Table 1, the electron cooling has to be used to reduce the normalized transverse 
emittance of the lower energy proton beam to 57rmm-mrad. Also, in that case 
the proton beam should have collisions only with the electron beam. Proton- 
proton collisions in the other two interaction points have to be avoided to allow 
for a higher proton beam-beam parameter. The maximum luminosity achieved 
in the present design is 4.4 x 10^^ cm~^s~^ in the high energy collision mode (10 
GeV leptons on 250 GeV protons). Possible paths to luminosities as high as 
10^^ cm~^s~^ are being explored, with studies planned to investigate the feasibil- 
ity of higher electron beam intensity operation. To achieve and maintain the Au 
normalized transverse beam emittances shown in Table 2, electron cooling of the 
Au beam will be used. For the lower energy tune of electron-Gold collisions, the 
intensity of the Gold beam is considerably reduced because of the reduced value 
of the beam-beam parameter limit for the electron beam. 

4-2 eRHIC: linac-ring design 

A linac-ring design for eRHIC is also under active consideration. This configu- 
ration uses a fresh electron beam bunch for each collision and so the tune shift 
limit on the electron beam is removed. This provides the important possibility 
to attain significantly higher luminosity (up to 10^^cm~^s~^) than the ring-ring 
design. A second advantage of the linac-ring design is the ability to reverse the 
electron spin polarization on each bunch. A disadvantage of the linac-ring design 
is the inability to deliver polarized positrons. The realization of the linac beam 
is technically challenging and the polarized electron source requirements are well 
beyond present capabilities (259). 
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Table 1: Luminosities and main beam parameters for e -p colhsions 
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Table 2: Luminosities and main beam parameters for e -Au collisions 



Figure 30 shows a schematic layout of a possible linac-ring eRHIC design. A 450 
mA polarized electron beam is accelerated in an Energy Recovery Linac (ERL). 
After colliding with the RHIC beam in as many as four interaction points, the 
electron beam is decelerated to an energy of a few MeV and dumped. The energy 
thus recovered is used for accelerating subsequent bunches to the energy of the 
experiment. 

4-3 Other lepton-ion collider designs: ELIC 

A very ambitious electron-ion collider design seeking to attain luminosities up 
to 10^^cm~^s^^ is underway at Jefferson Laboratory (256). This Electron Light 
Ion Collider (ELIC) design is based on use of polarized 5 to 7 GeV electrons in a 
superconducting ERL upgrade of the present CEBAF accelerator and a 30 to 150 
GeV ion storage ring (polarized p, d, ^He, Li and unpolarized nuclei up to Ar, all 
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Figure 30: Schematic layout of a possible linac-ring eRHIC design. 



totally stripped). The ultra- high luminosity is envisioned to be achievable with 
short ion bunches and crab-crossing at 1.5 GHz bunch collision rate in up to four 
interaction regions. The ELIC design also includes a recirculating electron ring 
that would help to reduce the linac and polarized source requirements compared 
to the linac-ring eRHIC design of section 4.2. 
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Figure 31: Schematic layout of a possible electron light ion collider at Jefferson 
Laboratory. 



The ELIC proposal is at an early stage of development. A number of technical 
challenges must be resolved, and several R&D projects have been started. These 
include development of a high average current polarized electron source with 
a high bunch charge, electron cooling of protons/ions, energy recovery at high 
current and high energy, and the design of an interaction region that supports 
the combination of high luminosity and high detector acceptance and resolution. 
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5 Detector Ideas for the EIC 

The experience gained at HERA with the HI and ZEUS detectors (260) pro- 
vides useful guidance for the conceptual design of a detector that will measure 
a complete event (47r-coverage) produced in collisions of energetic electrons with 
protons and ions, at different beam energies and polarizations. The HI and 
ZEUS detectors are general-purpose magnetic detectors with nearly hermetic 
calorimetric coverage. The differences between them are based on their approach 
to calorimetry. The HI detector collaboration emphasized the electron identifica- 
tion and energy resolution, whereas the ZEUS collaboration puts more emphasis 
on optimizing hadronic calorimetry. The differences in their physics philosophy 
were reflected in their overall design: HI had a liquid Argon calorimeter inside the 
large diameter magnet, whereas ZEUS chose to build a Uranium scintillator sam- 
pling calorimeter with equal response to electrons and hadrons. They put their 
tracking detectors inside a superconducting solenoid surrounded by calorimeters 
and muon chambers. HI placed their tracking chambers inside the calorimeter 
surrounded by their magnet. In addition, both collaborations placed their lu- 
minosity and electron detectors downstream in the direction of the proton and 
electron directions, respectively. Both collaborations added low angle forward 
proton spectrometers and neutron detectors in the proton beam direction. 

Both detectors have good angular coverage (approximately, 3° < ^? < 175°, 
where the angle is measured with respect to the incoming proton beam direction) 
for electromagnetic (EM) calorimetry, with energy resolutions of 1-3% and for 
electromagnetic showers a/E ~ 15% / ^/E~{GeV) + 1%. Hadronic energy scale 
uncertainties of 3% were achieved for both, with some differences in the a/E for 
hadronic showers, which were ~ R/\/E (GeV) + 2%, where R = 35% and 50% 
for ZEUS and HI, respectively. With the central tracking fields ^ 1.5 T covering 
a region similar to the calorimetric angular acceptance, momentum resolution 
a/pT < 0.01 (GeV) was generally achieved for almost all acceptances, except 
for the forward and backward directions. These directions were regarded at the 
beginning as being less interesting. However, the unexpected physics of low x 
and low (including diffraction in e-p scattering) came from this rather poorly 
instrumented region. And since the luminosity upgrade program, the low /?* 
magnets installed close to the interaction point to enhance the luminosity of e- 
p collisions (HERA-II) have further deteriorated the acceptance of detectors in 
these specific geometric regions. 

The EIC detector design ideas are already being guided by the lessons learned 
from the triumphs and tribulations of the HERA experience. All advantages 
of the HERA detectors such as the almost 47r coverage and the functionality 
with respect to spatial orientation will be preserved. The EIC detector will have 
enhanced capability in the very forward and backward directions to measure con- 
tinuously the low X and low regions that are not comprehensively accessible 
at HERA. The detector design directly impacts the interaction region design and 
hence the accelerator parameters for the two beam elements: the effective interac- 
tion luminosity and the effective polarization of the two beams at the interaction 
point. Close interaction between the detector design and the IR design is hence 
needed in the very early stage of the project, which has already started (255). It 
is expected that the detector design and the IR design will evolve over the next 
few years. The c-p and e-A collisions at EIC will produce very asymmetric event 
topologies, not unlike HERA events. These asymmetries, properly exploited. 
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allow precise measurements of energy and color flow in collisions of large and 
small-a; partons. They also allow observation of interactions of electrons with 
photons that are coherently emitted by the relativistic heavy ions. The detector 
for EIC must detect: the scattered electrons, the quark fragmentation products 
and the centrally produced hadrons. It will be the first collider detector to mea- 
sure the fragmentation region of the proton or the nucleus, a domain not covered 
effectively at HERA. The detector design, in addition, should pose no difficulties 
for important measurements such as precision beam polarization (electron as well 
as hadron beam) and collision luminosity. 

The EIC detector design will allow measurements of partons from hard pro- 
cesses in the region around 90° scattering angle with respect to the beam pipes. 
This central region could have a jet tracker with an EM calorimeter backed by 
an instrumented iron yoke. Electrons from DIS are also emitted into this region 
and will utilize the tracking and the EM calorimctry. Electrons from photo- 
production and from DIS at intermediate and low momentum transfer will have 
to be detected by specialized backward detectors. With these guiding ideas, 
one could imagine that the EIC barrel might have a time projection chamber 
(TPC) backed by an EM Calorimeter inside a superconducting coil. One could 
use Spaghetti Calorimetry (SPACAL) for endcaps and GEM-type micro-vertex 
detector to complement the tracking capacity of the TPC in the central as well as 
forward/backward (endcap) regions. This type of central and end-cap detector 
geometry is now fairly standard. Details of the design could be finalized in the 
next few years using the state of the art technology and experience from more 
recent detectors such as BaBar at SLAC (USA) and BcUc at KEK (Japan). To 
accommodate tracking and particle ID requirements for the different center-of- 
mass energy running (y/s = 30 — 100 GeV) resulting at different beam energies, 
the central spectrometer magnet will have multiple field strength operation ca- 
pabilities, including radial dependence of field strengths. Possible spectrometry 
based on dipole and toroidal fields is also being considered at this time. 

The forward and backward regions (hadron and electron beam directions) in 
e-p collisions were instrumented at HERA up to a pseudo-rapidity of 3. A 
specialized detector added later extended this range with difficulty to r/ ~ 4. 
Although acceptance enhancement in the regions beyond r/ = 4 is possible with 
conventional ideas such as forward calorimetry and tracking using beam elements 
and silicon strip based Roman Pot Detectors (181), it is imperative for the EIC 
that this region be well instrumented. A recent detector design for eRHIC de- 
veloped by the experimental group at the Max-Planck Institute, Munich, accom- 
plishes just this (136) by allowing continuous access to physics up to 6. The 
main difference with respect to a conventional collider detector is a dipole field, 
rather than a solenoid, that separates the low energy scattered electron from the 
beam. High precision silicon tracking stations capable of achieving Ap/p ~ 2%, 
EM calorimetry with energy resolution better than 20%/\/^, an excellent c/tt 
separation over a large range, all in the backward region (in the electron beam 
direction) are attainable. In the forward region, the dipole field allows excellent 
tracking and a combination of EM and hadronic calorimetry with 20% / ^/E(GeV) 
and 50%/y^E{GeV) energy resolution, respectively. This allows access to very 
high X ~ 0.9 with excellent accuracy. This region of high x is largely unexplored 
both in polarized and in unpolarized DIS. A significant distance away from the 
EIC central detector and IR, there may be Roman Pots, high rigidity spectrome- 
ters including EM calorimetry and forward electron taggers, all placed to improve 
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the measurement of low angle scattering at high energy. 

Although significant effort will be made to avoid design conflicts, the conven- 
tional detector using a solenoid magnet and the one described above may not 
coexist in certain scenarios being considered for the accelerator designs of the 
EIC at BNL. The main design line, presently the ring-ring design, may be par- 
ticularly difficult with only one IR. Options such as time sharing between two 
detectors at the same IR with the two detectors residing on parallel rails may 
be considered. In the case of the linac-ring scenario, several other options are 
available. Because the physics of low x and low does not require a large 
luminosity, nor is presently the beam polarization a crucial requirement for the 
physics (136). an interaction point with sufficient beam luminosity would be pos- 
sible with innovative layouts of the accelerator complex. These and other details 
will be worked out in the next several years. Depending on the interest shown 
by the experimental community, accelerator designs that incorporate up to four 
collision points (while still allowing two hadron-hadron collision points at RHIC) 
will be considered and developed. 
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